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ABSTRACT 


Atomic spectrochemical data is most commonly 
measured with a dispersive instrument. The dominant 
types of dispersive instrumentation are: the mono- 
chromator photomultiplier tube combination, the spec- 
trograph and the direct reader. The desired charact- 
eristics of a simultaneous multielement measurement 
system are present in the above types, but not one 
single type has all of them. The desirable charact- 
eristics are: fast, linear,wide dynamic range response; 
high sensitivity; wide simultaneous spectral coverage 
with good resolution; and direct conversion of the light 
intensity to an electronic signal. 

Attempts to achieve these characteristics have 
been based mainly on variations of the dispersive theme. 
DecOtablyrdriErerent approach is to multiplex the 
information with a Michelson interferometer. 

This is the Fourier transform spectroscopy 
technique which has been very successful in the mid and 
far-IR. To extend its use to the shorter wavelengths of 
the UV-visible region, the sampling of the interferogram 
must be considered because the clock signal (derived from 
a He-Ne laser) does not provide a high enough rate to 
prevent aliasing. The effects of apodization are examined 
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and a 2,5312 um sampling interval, 
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interpolated using the zero-filling 
technique. 
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Background and polystyrene spectra 
in the mid-IR. 
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CHAPTER I 
SPECTROCHEMICAL MEASUREMENTS 


The measurement of atomic spectrochemical data is 
one of the major techniques for acquiring qualitative and 
quantitative information about a sample. In an endeavour 
to obtain this information as accurately and as quickly as 
possible, there has been a continuing long range effort in 
many laboratories towards the development of effective 
simultaneous multielement analyses. A key aspect’ in this 
development is the design of spectrochemical measurement 
systems capable of simultaneously measuring spectral 
information over a wide range of wavelengths. This study 
1s)directed tosthat goalliby Anvestigating the feasibility 
of a Michelson interferometer Fourier transform spectroscopy 
approach, 

Some background is needed to put this technique in 
perspective. The spectrochemical measurement process is 
shown in general terms in Figure 1. Spectrochemical 
information about a sample is encoded by various processes 
such as emission, absorption, fluorescence, scattering and 
reflection into electromagnetic radiation with the information 
represented by the parameters of frequency, intensity, time, 
and space. This must be decoded to get useful analytical 


data about the sample. 


in af 2YEb feodimmitepacsoane aha :* Spates Sxl wel - 
ry 7 , 
ae 
ERA era 6 int wpignae Po err ae el eens i) ied Sa) sok ee: 


: ; Lone 
ti wiaokin, 2 br “hath ice sb (ee sole an Pie gS 


fre ovitns {tds et stb se to¥. haopeede= verifooil, att sin 


’ a. 
ny A iis 


a 


Hi PGs Sacay 4A aL, apr enearies eee peat STD ae ss | 
ditties Wlis Lo oderaarolerah tht | Chewes Be 

aT be Mls aeRO A! ESET AS ORS ah Shahn eae 
yaieaee cn Latinsiserssegs 3 ais nistcet 


teats gnisivessm. ys sujot aPALhe Sty tee bans ‘ae 

Chie aay verlignaleve e i 60 EO “Tae Pvt 8 
zitRtibeaad sit .ederet teed? “yO Jen eae ey beet aekb | 
WoGUeALISsTe GG ewes se) TALI TS2900T cae ee 162 fannie “a 


hae 


bedoss “ins Gud os bovsat ag oem gMoel _ 


~ 


al SepeCaTd PHStiStseovi Diseases ya pacty aay sae 
Initio: 41 Sau x rt wa i? Jaavoon 2k 

BeshGoore BUOL “EY Yd) cebosay Gi Lens & dpicin no) sadn 

) 

baa gaizetiace. .sdpgoesronls, wetaqyeeds. wiorecing af 
faistermoiens etd datw oolaskbat pitongaivaoels eagt eokioe 
ime Wefan, .Vonaipet} Bo eto wages, oly Yel baleen ot 

ea, beoddyl brs iiseg (392 09 bobaesh ad Jeu ein 


site only 


*wei8eIp YOOTq qUSWeANsesU TeoTWeYyooaqoeds °T eansTty 


NO!153143a S$ 
JDVdS V ONINSLLVOS 9 


viva IWIL © : 
IWOILAIVNY ‘ ALISNILNI Z echt : 
ADNINOG3IJ | 


NOISSIWS | 
pNouwwyosni], fF Nonviava |g __JNOWYWAOSNI 


LNOGV33 a= IVOILAIVNY (Qa ie 


__3G09N3 _ __ 


=~ el ces 
RO Piek4 ag 


By far the dominant technique for the measurement 
of spectra is the dispersive system based on the diffraction 
grating combined with a photomultiplier tube (PMT) or 
photographic plate detector. The main types of instrumen- 
tation for decoding the information, shown schematically in 
Figure 2, are the monochromator with a PMI, the spectrograph 
and the direct reader. 

The spectrograph has a photographic plate placed at 
the focal plane as the detector. Although the photographic 
plate is Capable Gf recording thousands: of lines in a single 
exposure aS a permanent record, it has a nonlinear intensity 
response, limited dynamic range and very tedious readout. 

In the monochromator PMT combination, the spectral 
information is dispersed along a focal plane where an exit 
slit allows one resolution element to reach the photomul- 
tiplier tube. The spectrum is obtained by scanning the 
dispersive element so that the frequencies reach the PMT 
sequentially. This can be very slow if high resolution 
and/or wide spectral coverage is needed. However, the 
PMT has the advantages of having very good linearity, wide 
dynamic range, high sensitivity and electronic output 
(light intensity is directly transduced into an electronic 
signal). 

The direct reader is an attempt to retain the advan- 
tages of the PMT and also achieve simultaneous spectro-~ 


chemical measurements, Here an exit slit-PMT combination 
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is mounted at each point in the exit focal plane of a 
dispersive instrument at which measurements are to be made, 
But even at best, only a small fraction of the spectral 
information available in the spectrum can be measured, 
Physical size and cost considerations limit the direct 
reader to about 60 channels, 

Over the last few years much effort has gone into 
the development of instrumentation that has the best 
characteristics of the above types. Desirable character- 
istics are: fast, linear, wide dynamic range response; 
high sensitivity; wide simultaneous spectral coverage with 
good resolution; and direct transduction of the light in- 
tensity to an electronic signal. 

The recent technological advances in solid state 
electronics have produced electronic image sensors as 
detectors (1,2) which over a moderate wavelength region 
provide a continuous multichannel measurement system. 

Such devices, which include silicon vidicons (3-6), 

secondary electron conduction vidicons G7) and=siiicon 
photodiode arrays (8,9), have recently been applied to a 
variety of multichannel spectrochemical measurements. 

At present these devices cannot match the photomultiplier tube 
performance in sensitivity and wide dynamic range. 

In addition to these multichannel systems, rapid 
scanning sequential systems have also proved useful Lor 


multielement analysis. Both a mechanical slew scan system 
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based on a programmable monochromator (10), and what amounts 
to an electronic slew scan system on an image dissector 
PMT (11) have been developed. 

Another potential approach to the overall problem 
is to dispense with the dispersive system completely and 
use a multiplex technique. Here the spectral information 
is encoded so that a single detector can be used to 
simultaneously measure a wide wavelength range (Figure 3). 
The most common examples of multiplex techniques are 
Fourier transform spectroscopy and Hadamard transform 
spectroscopy (12). In the latter technique, a dispersive 
element is still necessary but the spectral information 
is encoded in a binary code based on Hadamard matrices (13). 
In Fourier transform spectroscopy, the information is encoded 
in sine and cosine oscillations. This encoding is accom- 
plished by a Michelson interferometer. The encoded signal, 
which is measured with a single detector, is called an 
interferogram, It is necessary to take the Fourier 
transform of the interferogram to decode it and obtain 
the spectrum, Although Fourier transform spectroscopy 
has been associated with infrared measurements, there is 
no fundamental reason why it cannot be used for atomic 
spectrochemical measurements. 

An introduction to Fourier transform spectroscopy 
will discuss some advantages of the technique over conven- 


tional methods. Also, the basics of how the Michelson 
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interferometer operates and the method to obtain the final 
spectrum from the interferometer output will be explained. 
It is important to understand the additional factors which 
need to be considered when extending the use of Fourier 
spectroscopy to measuring atomic emission spectra in the 
UV-visible spectral region. The sampling of the inter- 
ferogram signal is an important consideration. Since the 
information in the UV-visible is at a higher spectral 
frequency than in the IR, a problem arises because the 
clock signal (derived from a He-Ne laser) does not 
sample the interferogram at a high enough rate to prevent 
aliasing. The output spectrum may not have the desired 
instrumental line shape. The process of apodization 
affects the line shape; therefore, apodization and its 
effects will be discussed. A major problem in Fourier 
spectroscopy is the phase error in the interferogram. 
The standard procedures for phase correction in FT-NMR 
and IR do not work effectively for line emission spectra. 
This problem and its solution will be discussed. 

A Michelson interferometer system has been built 
to study its use as a spectrometer. A block diagram is 
shown in Figure 4. The philosophy behind the design was 
to make a general purpose spectrometer that can be used in 
a variety of instrumental setups. Thus, the Michelson 
interferometer Fourier transform spectrometer can be used 


to take measurements in the ultraviolet, visible or 
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infrared regions, A variety of sources and detectors can 
be accomodated to suit the measurement. The mechanical 
design of the interferometer will be described along with 
the control circuitry and the measurement electronics. 
The software necessary to acquire and process the inter- 
ferogram will be described, A complete listing of the 
computer programs along with a detailed description is 
given in the Appendix. 

Finally some measurements taken with the Fourier 
transform spectrometer will be given to show the perfor- 
mance of the system. Spectra in the visible region were 
obtained with a 1P28 photomultiplier tube of hollow 
cathode lamps of several elements including Cr, Al, Mg, 
Ca, and a multielement lamp of Fe, Cr, Ni, Mn, Co, Cu. 
Flame emission spectra were also measured for Li, Na, Ba, 
Sr, and Ca, The Mg spectrum was measured in the UV with 
a solar blind PMT. Near-IR measurements were taken with 
a silicon photodiode detector of the flame emission of the 
aikali metals: Di, Na, K, Rb, Cse°s Quantitative aspectsof 
the system were studied with Li flame emission. Mid-IR 
measurements were taken with a TGS detector and a globar 


source and polystyrene film sample. 
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CHAPTER ITI 
FOURIER TRANSFORM SPECTROSCOPY 


A. Reasons for Using Fourier Transform Spectroscopy 


The question immediately arises as to what advan- 
tages the Fourier transform approach offers over the 
standard instrumentation and then what are some of its 
disadvantages, 

Fourier transform spectroscopy has been extensively 
used for spectral measurements in the infrared regions, but 
so far it has found little application for atomic spectro- 
chemical measurements in the UV-visible spectral regions. 
Much of the incentive for the development of infrared 
Fourier spectroscopy came from the promise of achieving 
two advantages brought to attention by P. Jacquinot and 
P, Fellgett: the throughput and the multiplex advantages 
which carry, respectively, their names in honour of their 
contributions to this field. These advantages result in 
the Fourier transform spectrometer making more efficient 
use of the light available from the source. 

The Jacquinot or throughput advantage (14) for the 
interferometer arises because the energy throughput, I, 
for a given resolution, R, is greater for an interferometer 
than a grating spectrometer. The reason is the 


larger acceptance angle for the interferometer without 
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degrading the resolution. The throughput for a grating 
instrument is limited by the necessity for an entrance slit. 


The derivation (15) gives the advantage as: 


ey T interferometer mane! Gaby. 
BRR ie a Leon duntea come ran ati turd nay ECM IR IE 
° “grating ) 


where l-equals the slit height and £ the focal length 
of the collimator of the grating instrument. A typical 
value for (1/f) is 1/50 (slit height of 1 cm and a 0.5 m 
collimating focal length). Thus the interferometer, in 
theory, has over two orders of magnitude nore energy 
throughput. The Jacquinot advantage is applicable to both 
infrared and the UV-visible regions. 

The multiplex (Fellgett's) advantage arises from 
the fact that the interferometer is observing all the 
spectral elements simultaneously whereas the scanning 
monochromator observes each spectral element sequentially. 
Suppose the eee nealigeet on of interest is between Ma and 
V5 with a resolution of 6v. The number of resolution 
elements M is then: 

M = (v, - v,)/8v ‘ (29) 

In a grating instrument each resolution element $v is 
observed for a time T/M, where T is the required time to 
scan between v, and Voe The signal received is proportional 
to T/M. Assuming that the noise is random and independent 
of the signal intensity, the signal noise is proportional 


to (t/a) t/2, Thus, the signal-to-noise (S/N) ratio 
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would be: 
S/N oc. (t/a) }/2 (3) 
For the interferometer each resolution element is observed 
for a time T SO that the signal received is proportional to 
T/2 (the factor of 2 arises because in a Michelson inter- 
ferometer half the light is reflected back to the source), 
The signal noise, assuming that noise is random and signal 
independent, is proportional to ql/2, Thus, for an 
interferometer the 
s/n oc TH/2 79 ; : (4) 

Comparing equations (3) and (4), the S/N for the inter- 
ferometer is better than for the monochromator by mi/2 79, 

This is the case in the infrared enon where the 
detectors are the major source of noise, and thas the 
noise is signal independent. However, in the UV-visible 
region the detectors are much more sensitive and the noise 
is determined by the fluctuations in the number of photons 
arriving at the detector, There can be several sources of 
this fluctuation in the number of photons (16,17). One 
source is the statistical nature of light. The noise level 
due to this is proportional to the square root of the 
source intensity. An even larger source of noise, especially 
when the source is a flame, may be source flicker, Noise 
due to source flicker is directly proportional to the source 
intensity so that this may be the dominant noise source 


at high signal intensities, Therefore, the benefit of 
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increased signal level for a resolution element due to the 
observation for the entire time T is cancelled because the 
noise level now increases as well. Hence, in the UV-visible 
region, the multiplex advantage is not realized. 

But the situation is not quite so clear cut. Both 
the nature of the spectrum and the distribution of the 
noise in the spectrum can affect the existence of a 
multiplex advantage in signal noise limited measurement 
situations, There can be a multiplex advantage with spectra 
of low density, such as atomic spectra. The advantage is 
realized if there is little or no broadband background 
radiation present or if the spectrum does not contain 
large numbers of much more powerful lines. This has been 
discussed by Filler (18), Hirschfeld (19), and Winefordner 
(20) and essentially the same conclusion has been reached 


for Hadamard transform spectroscopy (21). 


Another factor that must be considered is the possi- 
bility of a so-called multiplex disadvantage (piles ose ee a SS: 
arised primarily because the distribution of photon noise 
in a spectrum measured with a scanning instrument is dif- 
ferent from that in one measured with a Fourier spectrometer. 
In the scanning case, the rms (root mean square) noise is 
greatest where the signal is greatest CEhateis, eat Che top 
of the spectral peaks). In the Fourier case, the noise tends 
to be spread out more or less uniformly throughout the 


entire spectrum. Thus, the S/ nN ratio of strong peaks 
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should improve, but weak spectral lines may be obscured by 
the noise from strong lines that ends up distributed along 
the baseline of the spectrum. Therefore, if the spectral 
peaks of interest are small compared to peaks that are not 
of interest or background radiation, there is a decrease in 
the S/N ratio compared to the scanning instrument. This 
problem also exists in Hadamard transform spectroscopy (21), 
and some experimental verification of this multiplex dis- 
advantage has been reported (22). 

Thus, at this time no definitive conclusions can 
be reached about the existence and importance of 
multiplex advantage and/or disadvantages in photon noise 
limited multiplex spectrochemical measurements. More 
experimental work is necessary to clear up the situation. 
What the above discussion does mean, however, is that the 
promise of any substantial multiplex advantage cannot be 
a driving force in extending Fourier transform techniques 
into the UV-visible region. In fact, some disadvantages 
may exist with respect to S/N in certain measurement 
situations. | 

However, signal-to-noise ratio is not always the 
only and overriding consideration when carrying out a 
spectrochemical measurement, It becomes important only 
when nearing the limits of detection and many measurements 
are made well above this limit. The impressive success 


and capabilities of the Fourier transform technique in 
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the infrared have shown that some important advantages 
result from the nature of the instrumentation used to 
implement Fourier transform spectroscopy. These 

advantages are not dependent on the existence or realization 
of the multiplex advantages. 

Some of these advantages are: spectra can be 
measured with a very accurate and precise wavenumber axis 
which is predetermined by the sampling interval (no 
calibration necessary), a wide range of frequencies can 
be simultaneously covered with proper utilization of 
aliasing, high resolution can be achieved in a relatively 
compact system along with speed of operation, the 
resolution function is easily controlled and manipulated 
as an inherent step in data reduction by use of apodization 
techniques, and computerization of the spectrometer is 
facilitated, 

The digitization of the interferogram in a Fourier 
transform spectrometer is normally controlled with a 
clock signal derived from a He-Ne reference laser. This 
provides a final Spectrum which has a very accurate 
wavenumber axis with the values along the axis easily 
calculated from the wavenumber of the laser line (no 
external calibration needed). This inherent accuracy 
of the axis greatly facilitates intercomparison of 
digitized spectra for small peak shifts and peak shape 


perturbations. Also, the accurate axis makes 
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identification of spectral peaks easier. 

The spectra are in a digitized form as an inherent 
part of the data reduction, that is, the spectra are stored 
in the computer as part of the data processing. Hence, 
further processing such as ratio, peak subtraction, and 
correlation of spectra is easily accomplished by software. 
A benefit of having the interferogram stored in a 
computer is that, if the resultant spectrum does not have 
the desired resolution function, it is a simple matter to 
process the interferogram again with a different apod- 
ization function (23). Another possibility is to apply 
some form of digital filtering to smooth the spectrum or 
enhance resolution (24,25). 

A Fourier transform spectrometer is capable of very 
high resolution (26-28) while maintaining relatively high 
throughput and with a relatively compact system. A 
grating instrument for high resolution can be quite 
bulky. For example, a 2 meter spectrograph is rather 
large and heavy. On the other hand, an interferometer 
with 8 cm mirror movement (16 cm optical retardation) can 
achieve nominal resolution of 0,0625 oa (0.00225 nm at 
600 nm) and still easily be placed on a tabletop. 

With proper aliasing, a Fourier transform spec- 
trometer can simultaneously cover a wide range of wave- 
lengths while maintaining relatively high resolution. 


This aliasing advantage will be considered in Chapter III. 
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Considered in the context of capability, both the 
cost and basic simplicity of Fourier transform spectroscopy 
can be considered advantages. Present commercial 
instruments have not supported the cost advantage since 
they are very expensive, but without question, effective 
systems can be made at considerably lower costs. The 
simplicity is present in the fact that there is but one 
moving part in the spectrometer, the moving interferometer 
mirror. 

As a result of these additional considerations, it 
seems that a Fourier transform approach is feasible in the 
UV-visible region (29,30) as well as in the infrared where 
it has without doubt shown very fine performance. The 
types of instrumentation and many applications of the 
Fourier transform approach has been discussed in the 
literature (31-37). One major area where Fourier 
transform spectroscopy has been used is in astronomy where 
much of the development for high resolution work has taken 
place (26,27,38-42). The basis of the Fourier transform 


technique is summarized in the next section. 
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B. The Michelson Interferometer 


The main optical component of most Fourier trans- 
form spectrometers is the Michelson interferometer, A 
simplified diagram is shown in Figure 5, The inter- 
ferometer encodes the spectral information from the source 
via two beam interference, The output signal at the 
detector is called an interferogram. 

The most common arrangement of the Michelson inter- 
ferometer is with a beamsplitter at 45° to the source. 
In the visible and infrared regions the beamsplitter 
material is too thin to be self-supporting, thus, it is 
coated onto a support material. A second equal thickness 
of the support material, called the compensator, is placed 
in one arm of the interferometer to equalize the optical 
path length in both arms, The beamsplitter ideally 
reflects 50% of the light and transmits the remaining 50%. 
The two divided beams of light reflect off plane mirrors 
at right angles to one another and are recombined by the 
beamsplitter, Half of the recombined light goes to the 
detector and the other half is returned to the source. 
One of the plane mirrors is fixed and the other mirror 
can be moved in a direction perpendicular to the mirror 
surface, 

To see how the interferometer encodes the incident 


radiation, assume for simplicity that there is monochromatic 


in 


‘] 


ome: 


edie tag ieta) ee 2a? 


'oAskh eid’ bi of sashimi 
") U i i'n ° 
Pri? yi 6 ma eu ara is Pp ie. 2 . 
i) feteo@e0 er carat Saae 
r' i 5 meas pi rues Ne: 
») 
re eh i sade 
- 
2 r 
- vo 
. aT i a ) oes ‘ce 
i te iy | mie 
: ; 
ats 


SOURCE 


BEAMSPLITTER 
COMPENSATOR 
—Be DETECTOR 


MOVABLE 
MIRROR 


FIXED MIRROR 


Figure 5, Schematic diagram of a Michelson 
interferometer, 
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light from the source. When the optical path lengths 
of the two arms are equal, the two beams of 
light are in phase when they return to the beamsplitter 
and thus interfere constructively. If the movable mirror 
is offset z of the wavelength of the incident light, the 
optical path lengths will differ by % a wavelength. When 
the two beams are recombined at the beamsplitter, they are 
now 180° out of phase and will interfere destructively. 
As the mirror is continuously moved, the recombined beams 
of light will alternately interfere constructively and 
destructively, resulting in a cosinusoidally changing 
output at the detector. The maxima of the output 
(constructive interference) occurs for every % ae as) 
wavelength movement in the movable mirror, or, in other 
words, for every one wavelength difference in optical 
path length of the two arms. An example of a monochromatic 
source is a He-Ne laser and its interferometer output is 
shown in Figure 36(a) (lower trace). 

The frequency of the output cosine wave depends on 
the wavelength of the incident radiation and the velocity 
of the movable mirror. Thus a very high optical frequency 


s Hz) is uniquely encoded in 


(for example 600 nm is 5 x 10 
the form of a low frequency oscillation. (In ithe case (of 
the interferometer used in this study, the 632.8 nm He-Ne 
laser line is modulated to around 4 kHz. 


If the incident radiation is broadband, there are 
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many wavelengths of light. Each wavelength is modulated 
to a different frequency. The output of the detector is 
then the summation of all the frequencies present. The 
output with a white light source (tungsten bulb) is shown 
in Figure 6. At zero path difference, ZPD, (optical path 
lengths equal in the two arms of the interferometer) all 
the frequencies are in phase giving a large output signal, 
As the mirror is moved away from the zero position, the 
frequencies rapidly sum out to a steady average value. 

In Figure 6 the white light interferogram was digitized 
at optical path difference intervals equal to one wave- 
length of the He-Ne laser 632.8 nm line, The ZPD 
position occurs at about point 2048 since the 4096 point 


interferogram was sampled on both sides of the ZPD position. 
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Figure 6, Example of an interferogram of a 


broadband source (tungsten bulb). 
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C. The Fourier Transformation 


The basic process for decoding the interferometer 
output is to convert the amplitude--path difference 
information represented in the interferogram to the 
amplitude--frequency domain of a spectrum. This process 
is called the Fourier transformation. 

The detailed theory and mathematics of the Fourier 
transform has been extensively documented in the literature 
(15,31, 32,34,43-48) and will not be reproduced here except 
in summary. 

Mathematically the ac signal called the interferogram 


can be expressed as 


+ 0o 


Gg) J B(v) cos(2txv) dv (5) 
-00 
where I(x) is the intensity of the output as a function of 
mirror movement x; B(v) is the intensity of the source as 
@pLunctlOnsOfewavenumbersy, jehatwis uber spectrum : 
Equation (5) is one-half of a cosine Fourier transform 
pair (45). The other half is 
+20 
By) = ii LCx) «cos(20xwWedx i (6) 
~ 00 
Equations (5) and (6) provides the relationship between 


the interferogram, I(x), and the spectrum, B(v). The 
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spectrum can be obtained by taking the Fourier transform 
of the interferogram, 

If all the optical components were ideal, the 
interferogram in Figure 6 would be perfectly symmetrical. 
This is not the case, therefore, there are sine components 
in the interferogram, Thus, in general, the interferogram 
and spectrum must be related by the complex Fourier 


transform pair, 


+00 


Vx) = | B(v) exp(2tixv) dv G7) 
BGA = I(x) exp(€2tixv) dx ‘ (8) 


The exponential term contains both cosine and sine 
contributions (45). 

In practice, the interferogram cannot be measured 
to infinite path difference, but only to a maximum path 
difference Rr ete Also, the interferogram is measured at 


Certain intervals of position; that is, the function I(x) 


is not continuous. The integral is then replaced by a 


summation. Equation (8) is then replaced by 


N 
Beye oe Len exp (2nivkh) (9) 


where h is the sampling interval and N=x_,, /h, N being 


the number of interferogram sample points from Xp to X%,.,- 
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The calculation of the spectrum by direct 
application of equation (9) is the classical method, 
However, it requires long computation times even by 
large computers, especially as N increases, because the 
time needed is proportional to Ne But since the 
development of the fast Fourier transform (FFT) 
algorithm by Cooley and Tukey (49), the most common 
method to calculate the spectrum is with the FFT routine. 
A very good description of the algorithm is given by J. 
Connes (50) and by others (44,51). 

In this study the FFT was considered as a "black 
box" to convert from the time to frequency domains. The 
actual routine used is the same as used by Horlick (52) 
and is written in floating point FORTRAN. The size of the 
input data array must be an integral power of 2, for 


eS OO CLL E PR Oaree nde balichtiok an 


example, 272512, 2 
integral power of 2, it can be filled with zeroes to the 
next power of 2. It should be mentioned that the speed 
of the program allows a 4096 point FFT to be performed in 
about two minutes on the PDP-11/10 minicomputer. A 
single precision machine language FFT routine available 
for this computer requires less than 10 seconds to 
transform the same number of data points. 

In the routine used, the odd numbered points of the 


input array, length N, are placed in the real input array 


and the even numbered points are placed in the imaginary 
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input array. The real “and imaginary input arrays are also 
of length N; therefore, since only the first N/2 points 
have been filled so far, the remainder of the arrays are 
Filled with zeroes. This produces a smoother spectral 

plot because two output points are generated per resolution 
Clemente Ot COUrSe 4. tNe spectral “resolution, «whichis 
determined by the maximum optical retardation in the inter- 
ferogram and the apodization (discussed in Chapter IV), is 
not altered by filling the arrays with zeroes. 

The output of the FFT routine consists of two series, 
X(J) and Y(J), which are the real and imaginary intensity 
components and J is the index of the frequency components 
(from 1 to N). The amplitude spectrum of these frequencies, 
tj us Calculated by taking the root sum of squares, Of 
these two series: 

ame () 1) eae (10) 
Normally only the amplitude spectrum is calculated, The 
phase spectrum can be calculated from the real and imag- 
inary series (23). 

The wavenumber of a peak in a spectrum is calculated 
using the J index corresponding to the point having the 
largest value in that peak: The caleulatiton is based on 
the He-Ne laser wavelength of 632.8 nm. The, precision in 
the peak location is shown experimentally to be within one 


Holmes (lablomy2 4) or 0.03. nmin ithesvacinitv. of 400 nm. 
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CHAPTER III 
SAMPLING THE INTERFEROGRAM AND THE 


CONSEQUENCES OF ALIASING 


In order to perform the digital Fourier transform 
of the interferogram, the analog signal must be sampled 
and digitized at appropriate intervals of time. The 
sampled interferogram should be an accurate representation 
of the continuous analog interferogram; that is, no 
information should be lost. The importance of accurate 
digitizing is obvious, A data point must be acquired at 
equal intervals of optical retardation. A missed or 
extra data point can have adverse effects on the spectrum 
(53). Also, the analog-to-digital converter must be accu- 
rate and have enough resolution. At least a 12-bit 
converter is necessary. For broadband type signals, even 
more bit resolution may be needed to accurately digitize 
the small fringes at large optical retardation. 

Not so obvious is the importance of the sampling 
interval. Information theory gives the answer to what 
should be the appropriate sampling time interval (45). 

The frequency content of an analog signal is zero above 


a certain frequency, f To sample the analog signal 


max ° 
properly, the time interval between samples must not 
exceed 1/2f,... Another way of saying this is that the 
sampling frequency must be two times the highest 


frequency present in the signal. If the sampling rate 
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is too low, high frequencies will be improperly sampled 
and will show up as spurious low frequencies. This is 
termed aliasing. At near-IR and shorter wavelengths, 
aliasing is an important consideration; therefore a 
thorough understanding is necessary to interpret the 
spectrum. A simple illustration may explain aliasing more 
clearly. 

In Figure 7(a), spectrum A resulted from the 
Fourier transform of a synthetic interferogram which 
contained three frequencies, each of different amplitude. 
The interferogram was sampled at a rate such that all 
frequencies, up to X were properly sampled. For example, 
if X=200 Hz, then the sampling rate must be at least 
400 Hz (a sampling interval of 2.5 msec). Then every 
other sampling point was dropped, the spectrum was cal- 
culated again, resulting in spectrum B, Since the 
sampling interval is now twice as long, 5 msec, only 
those frequencies up to Y (100 Hz) are properly sampled. 
Peaks 2 and 3 now appear as spurious low frequencies and 
their position can be predicted as a folding over about 
the central point of the original frequency axis. The 
175 Hz signal (peak 3) appears as a signal of 25 Hz; 
that is, it has an alias of 25 Hz. The 125 Hz signal 
(peak 2) has an alias of 75 Hz. Peak 1 remains at its 
proper position because it is still properly sampled. 


Figure 7(b) illustrates how a frequency becomes aliased. 
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Figure 7. Illustration of aliasing. (a) Aliasing of spectral 
peaks due to undersampling of the interferogram. 


(b) Aliasing of a 175 Hz sine wave to a 25 Hz 
sine wave by undersampling. 
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The 175 Hz sine wave sampled at a 400 Hz rate is shown 
labelled as A. Every other sampled point is omitted, 
simulating the undersampling in spectrum B, and the 
remaining points are connected (shown labelled as B), 
A 25 Hz sine wave is the result. A similar exercise 
with a 125 Hz sine wave would result in an alias of a 
75 Hz sine wave. 

It is important to note that if peak 1 were not 
present, the fold over of peaks 2 and 3 would not be 
serious as it occurs in an exactly predictable manner. 

Of course, the approximate frequency range of the signals 
must be known to predict whether the spectrum is aliased. 
Therefore, in practice, it is the bandwidth of the 

Signal and not just the highest frequency which 
determines whether a certain sampling rate is acceptable. 

In the previous example (Figure 7), if peak l 
did not exist and the bandwidth was known to be only 
100 Hz (from Y=100 Hz to X=200 Hz), a sampling rate 
of 200 Hz would still aliase the spectrum (as in spectrum 
B). However, the spectrum can be interpreted because an 
a priori knowledge of the bandwidth eliminates the 
possibility that peaks 2 and 3 are in the 0 to 100 Hz 
region. There is no overlap of regions so that a sampling 
rate twice the signal bandwidth is acceptable. But, 
care must be exercised in the proper labeling of the 


frequency axis. 
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Normally, in Fourier transform Spectroscopy, 
aliasing is avoided as the undersampled modulation 
frequencies show up as spurious spectral information 
(fold over). This can be very confusing when there is 
much broadband energy present. However, with line 
spectra, as commonly measured in atomic spectroscopy, it 
is possible to use aliasing to advantage. To be useful 
one must have an a priori knowledge of the sampling rate 
and the spectral region of the signal. The bandwidth of 
the signal must be known, but it can be determined from 
the instrumental components. Factors limiting the 
bandwidth are: optical transmission of the beamsplitter, 
beamsplitter support and compensator materials; 
reflectivity of the mirrors and/or lenses used in the 
system; source output; detector response, both spectral 
sensitivity and frequency response; detector amplifier and 
electronic filter bandwidth; and any optical filters 
placed in the light path to intentionally block out some 
spectral regions. 

A practical illustration of the use of aliasing is 
the measurement of the flame emission spectra of lithiun, 
potassium, rubidium and cesium with the Fourier transform 
spectrometer. Table I lists the emission lines of these 
elements in the near-IR region along with the region 
from Table II in which they fall. Table II will be 


explained shortly. Figure 8 gives a schematic diagram 
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TABLE I 


in the Near-IR 


Major Emission Lines of Na, Li, K, Rb, Cs 


Element 


Na 


Li 


Rb 


Cs 


eet 


589.00 
389.99 


670.78 


766.49 
769.90 


780.02 
794,76 


852.11 
894,35 


Region(in Table IT) 
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of a spectrum of these elements with the frequency axis in 
wavenumbers. The amplitude of the lines are arbitrary. 
In order to measure these lines of Li, K, Rb, and 
Cs. without aliasing, a bandwidth of 15,802,8 cm! is 
necessary. However, the standard He-Ne laser used as a 
reference has a wavelength of 0.6328 um (15,802.8 em7!), 
The clock signal for sampling is derived from the laser 
| fringes generated by the interferometer, A clock pulse 
is generated for every 0.6328 um of optical retardation; 
that is, the signal interferogram is sampled at 0.6328 um 
ie rales This means that the shortest wavelength of 
light that can be properly sampled without aliasing is 
1.2656 um (7901.4 aap To achieve the bandwidth of 
15,802.8 cm7!, the sampling interval must be 0.3164 um. 
With the 0.6328 um sampling interval, the bandwidth covered 
is 7901.4 cm-+, The 7901.4 to L5g8028 cm + region then 
becomes aliased and overlaps the 0 to 7901.4 em + region. 
This is shown in the bottom diagram of Figure 8. In this 
case the overlap is not serious since there is no signal 
in the 0-7901.4 om? region (due to the Si photodiode 
detector's insensitivity and the quartz beamsplitter not 
being efficient in the mid-IR region). Therefore, as 
long as one is aware of the aliasing, the frequency axis 
can be labelled properly and the spectrum can be 
interpreted. 


Now, what if there is spectral signal above 
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1,002.9 omy? An example is where sodium is present 
along with the previous four elements, shown schematically 
in Figure 9, Even if the sampling interval were 0.3164 um 
i605 oem) the Na doublec vat’ 599 nn would Stadt ve 
aliased., The overlap is shown in the middle diagram of 
Figure 9 where the Na doublet is aliased to a position 
between the K and Li lines. With the 0.6328 um sampling 
interval, the result is shown in the bottom diagram. 
The Na doublet has been aliased twice. There is spectral 
information from more than one aliased region and 
interpretation is less trivial. But, as long as the 
aliased lines do not overlap exactly, interpretation is 
possible. One possible method to determine which region 
a line belongs is to use an optical filter to block out 
the other region or regions and measure the spectrum again. 
The disappearance or continued presence of a line will 
locate the region in which that line belongs. Also some 
regions may be ruled out due to lack of detector response, 
Another example is diagrammed in Figure 10. The 
chromium triplets at 427 nm and 360 nm are shown 
unaliased in the top diagram. If the sampling interval 
were 0.3164 um, the triplets would be folded over into the 
0-15,802.8 emi region pivoting about the 15,802.8 se 
point (middle diagram). With a 0.6328 um sampling inter- 


val, the bandwidth is 7901.4 om + and the Cr 427 nm triplet 


-1 ; 
is folded over again about the 23,704.2 cm point. 
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This overlapping of spectral regions due to aliasing 
is tabulated in Table II based on the sampling rates 
derivable from the He-Ne laser 632.8 nm line. Column III 
is the overlapping regions for a 0.6328 um sampling 
interval which has a bandwidth of 7901.4 cm™/, This 
is termed the direct sampling rate since a data point 
is taken for every laser fringe (0.6328 um optical 
retardation). In order to avoid confusion, the notation 
used will be explained. The regions are numbered at the 
left of Table II. When reference is made, for example, to 
the 15,802.8-7901.4 em7! region under the direct sampling 
rate, that region will be called region (5-8). Not all 
the regions listed under one sampling rate will be 
present in an aliased spectrum, As mentioned before, 
some regions can be optically eliminated. Another 
method is by using low and high filters to electronically 
block out regions. For example, if the mirror drive 
velocity is set to give a laser fringe of 4 kHz, this means 
that an optical signal of 632.8 nm is modulated to a fringe 
frequency of 4 kHz and lower wavelengths (higher optical 
frequencies) are modulated to a higher fringe frequency. 

If the frequency response of the detector amplifier is 
limited by a low pass filter set at about 4 kHz, all 
optical signals which are modulated to a frequency greater 
than 4 kHz (all wavelengths less than 632.8 nm) will be 


eliminated. A high pass filter can also be set to 
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eliminate low optical frequencies, 

Clearly then, to work in the UV-visible region, 
either the sampling rate must be increased or aliasing must 
be tolerated. The frequency a the reference laser 
modulation can be increased by use of optical techniques 
such as double passage of the laser through the inter- 
ferometer or by electronic frequency multiplication using 
phase locked loops (54). However, the high sampling rate 
necessary to sample the UV-visible modulation frequencies 
without aliasing quickly results in a prohibitively large 
number of data points that must be digitized in the 
interferogram to achieve reasonable resolution. Thus, 
the alternative of aliasing the spectral information in 
the interferogram must be used. 

To show that electronic multiplication of the laser 
modulation frequency was possible, the reference clock was 
frequency quadrupled using phase locked loop techniques. 
The analog laser fringes were also used as the spectral 
signal, Figure 11 shows the spectrum from the Fourier 
transform of the sampled laser interferogram, The sampling 
interval was 0.1582 um (x4 sampling rate) because a 
data point was taken at every 1/4 wavelength of the laser. 
The bandwidth was 31,605.6 om7> (Column V in Table II) 
and the laser line should be exactly at the midpoint of 
the spectrum at 15,802.8 cm7*, The peak does occur at 


this position. Spectrum (a) is the output with no 


be 


ds a ; 


e eater i a ay ‘Stid 2 pytip iat 


causbtr ved Ls - rau 2 
antiad tents gous se wel ai Ye : sieves 


‘7 Aa eee pe eee ee 
nites ngiahatigts iu0 pone pea, oli soto 03 


ort ty ‘ani igi » (ps fee MA. mee, (aE) ‘equet 
eal pen #Gitphubom 440 tae eit? ic, at 8 
ay 58 Pn! he ‘i. 

Lenk sleet tktdowg es we Ra faeeD anaug grt wee 


SSABL Vey te tite 


ra 3 me Poy ee cere Pha le ee eoob ie 


? 


Hosa 4. en ape, iy ihied 


Sha £ % OL Tf! +h) 4 = etc rPaAnate at works ee ‘a 
- — > — _ é - & 
v a 

. A 1 * £ 3 

ga kooks sara7 bi et sas . aad" rere ‘onesie a 2 tjedt 


meBriy Une Sh) “ita 


a) ae - Ji Nie 
aagee Gi: Sia Non (Rast. Sea yes ale ns ad 
| shiek Quan mF 


- oa antici aP- ole e - al <A Lee ey. 7 
we sti! i: +t ae merits veh © ore Wa hae Waal wt Diet 
i rm nh a + soy on 
Snrie eo. Bris NE aa s . oO 
x. 
-% ‘ 
To & 4 ie ry < \ 56: : 2 
f "eet eS =o «2 7 ' 
‘SoG 4 ‘ fom Naor eel : ¢ VS & 
~ Ry > = 
(wy (« eT cf y pita: as) av ue . i 


forwabh ais Ja yisosx4 ad Hangs 2 


45. even Boob Seog af |. BP OM oe 


7 wih y. war dL 


42 


*uoTIezT pode 
te[Tnsuetaq (q) *uotjeztpode on (e) °*s3er Sut{[dues vy 
3e@ peTdwes wers0rejrzequT ‘feuTT um gz7c9°¢ APSeT eN-8H “TT vans ty 


(,4>) ADNINOIZY 
909 LE ZOe'Si | O- 


(9) 


Ce - ‘ 5 
Pe “i 
eh & 
7 
i 7 
ie p> ] 
‘a 
ti 
Y ’ i 
> ~ r | 
bee 1 
*) 
= | \ 
; Al 
| 
i 
a | 
- 
— 
i) 
“~h, 
} 
: 


apodization of the interferogram. Spectrum (b) is the 
result of triangular apodization before performing the 
Fourier transform. Apodization will be discussed in 
Chapter LV 

Electronic frequency doubling by the same technique 
gives a sampling interval of 0.3164 um. The interferogram 
of the flame emission of a Li, K, Rb, Cs solution was 
sampled at this rate (x2 sampling rate) and the resulting 
spectrum is shown in Figure 12(a). At this sampling 


1 


interval, the bandwidth is 15,802.8 cm™ and all the 


spectral lines are properly sampled. There is no aliasing. 
The direct sampling rate with a bandwidth of 7901.4 cm7+ 
results in spectrum (b) where all the lines have been 
aliased, These are the real spectra that were diagrammed 
in Figure 8, It should be noted that all the spectra 
measured with a Fourier transform spectrometer has an axis 
which is linear in frequency (wavenumber) and not linear 
in wavelength. This set of data was taken with the 
PDP=-8/e system which had a limit of a 512 point transform. 

With the limited number of data points and the 
wide spectral coverage (over 220 nm), the resolution is 
quite poor. The potassium doublet is not resolved at all 
and the rubidium doublet is only partially resolved in 
spectrum (a). In spectrum (b), the resolution is doubled. 
This is due to the fact that the sampling interval is 


twice as long; therefore, with the same number of data 
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K 
(a) 
Rb ' 
-t 
sie ve 
~s ; 
eS me sh 
ae | A— 7 8--> 
0) 7901 ern 15,803 cm 
K 
(b) is 
Li Mine 
| I\, Rb 
&s 
| | | Cs 
OLS |e eee A eae 2a 
r Flan perieaae  eat Fit El? ean eae 
15,803 cm"! 11,852.em! 7901 cm! 
K 
(c) 
Sea ae (Lithium) "7 (Potassium, Rubi hud 
15,803 cm’! 13,327 cm”! 11;852 cm 
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Figure 12, Spectra of Li, K, Rb, and Cs using 
a) X2 samplin rate, (b) direct, 
c) $2, and (d) #4 sampling rate. 
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points (512 in this case), the interferogram has been 
sampled to twice the optical retardation which determines 
the ultimate resolution, Table III tabulates the 
resolution that can be achieved with a triangularly 
apodized interferogram using different sampling intervals. 
The resolution is also given in wavelength units at 
770 and 280 nm. It should be noted that the resolution 
is constant in wavenumber for ane Situation; therefore, 
the resolution increases with decreasing wavelength. 
For example, 16 om7+ resolution at 770nm results in 
0.95 nm resolution; but, at 280 nm, the resolution is 
02125 -nm. 

If the sampling interval is 1.2656 um (a data 
point every other laser fringe), the bandwidth is 
BDU s ile The overlapping of the regions due to 
aliasing is tabulated in Column II of Table II. This 
is the $2 sampling rate. The spectrum is shown in 
Figure 12(c). The resolution is again doubled and now 
the potassium doublet is resolved. The cesium lines 
(in region 6) have aliased into the potassium-rubidium 
region (number 7). The region numbers are those given 
in Table II, 

If the sampling interval is made even coarser to 
2.5312 um, the bandwidth is 1975.35 em’. The regions 
overlap further (Column I of Table Tie) with all three 
regions overlapping via aliasing. The spectrum is shown 


in Figure 12(d). The theoretical resolution at the 
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TABEE EL I 


Resolution with a Triangularly Apodized Double-sided 


Interferogram of (a) 512 points, (b) 4096 points 


Sampling 


Rate cin 
X4 (a) 256 
X2 (a) 128 


Direct (a) 64 


Direct (b) 8 


Resolution 


mm at 770nm 


Oe. 
Tes) 
Beis 
way) 
CaS 
0.4/7 


nm at 280nm 


2.00 
1,00 
0.50 
0.25 
0) sAhas) 
0,063 
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different sampling rates is given in Table III for a 
triangularly apodized double-sided interferogram of 

512 points. Since the potassium doublet has a 3.4 nm 
separation, it is expected not to be resolved until the 
sampling interval was 1.2656 um. The spectra point 
this out. 

The flame emission was measured again when the 
new data system with the PDP-11/10 minicomputer was 
available, With 4096 point interferograms, the 
resolution is better with the direct sampling rate than 
with a #4 sampling rate and a 512 point interferogram. 
Figure 44 in Chapter VII shows a spectrum of the above 
four metals plus sodium from a 4096 point interferogram. 
The Na doublet has only a 0.59 nm separation, but it is 
resolved. Table III predicts this, pointing out that 
the resolution capability is fairly close to the calculated 
resolution, 

From the series of spectra in Figure 12, it is 
clear that aliasing can be used to advantage in 
optimizing spectral coverage and resolution when making 


Fourier transform atomic spectrochemical measurements. 
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CHAPTER IV 
APODIZATION AND ITS EFFECTS 


As mentioned earlier, the interferogram of a 
monochromatic source is a cosine wave. If the source 
were infinitely narrow, the cosine wave would continue 
£Or infinitely large optical retardation. However, in 
Praceice, this+is not observed. First, certain 
spectral or instrumental conditions result in the 
interferogram being damped as optical retardation 
increases. Factors contributing to this are natural 
line broadening and solid angle effects (15). Second, 
the interferogram is sampled to an optical retardation 
of Kee SndanoLeconinicanitye, his truncation ors tne 
interferogram is a type of self-apodization. The term 
apodization refers to the multiplication of the inter- 
ferogram by a relatively simple function. The effect of 
multiplying the interferogram by a function is equivalent 
to the convolution of the spectrum resulting from the 
Fourier transform of the interferogram with the Fourier 
transform of that function. This is shown in Figure 13 
where the asterisk "*" refers to convolution. The 
bottom diagrams describe apodization and the top diagrams 
describe convolution where the spectral peak is convolved 


with the instrumental line shape to give the output 


Spectrum, 


48 


cate a) cata 
a ; hs 


at wew Sith eson - ageoe, ts 


domiacds bittow SVB - aise aH¢ 


j 9, a 
ri 4 a WOH P35 eB pe | gis ae ee ve e ane bats vice ‘i 


. += ae ra) Weigelt " ‘gon at ats “ 


* ne - 2 > 
aT) < a . . 
. 
Ge 
_—— = ° - vy ay nS 
AAtisBhahiIst Sao 2 5 Perio. SiiLle 
ey ‘ R35 £ 


j . — a 7 4 
feunven ors sino od xis eOtiveee baat date 
So ee heh oe, et ~ , + 4 - ne " 

/ l t . 
byt a 7.) =e 7 a et ny a 
iGgvelc gi ' j ) : ; Lt oe 


‘be 


o™ t ae +¢ ‘ { =) 449) a! geet 5% £ 
73h F 7) ae hue Os: bah frtes of 
f F] ' 
’ a. — atin t » Oar bos dia | 
’ \ a) & a | k =a » - 
j 


1 i {pees wyysosqe See 
‘we 
aes _ &4 5) rl yy ‘BS _ oh ie "Gi { ‘ +") ry fy. STO She se a! 
| a. 


et surgi ni ywofi2 ied, |e’ ~ sro St Aei2 %e ooh 


wi 

sell! ; Pie . eyed 
< : = 
(ts @ be r . on 


eT a fie Loe) qas 3 rel vie fobissti os 


aed is hee ae tal % tn i Aa Vs 
bey lovncs; a 6 eae S12 S23Q- » = 


49 


“uoTqezrpode (7) °(9). *(p) *uor4anjoauoo (9) °*(q) *(2) 
sUOTIeZzTpode pue UCTINTOAUOS UseMJeq dTysuOoTAeTeEr sul “°ET san3Ty 


NOIHLDNAS ONIAFLOOWS 
LhelNO Wed G3zIIGOW NivWOG yalxno4 INdLNO 1V38 


NCOILVWYOSSNVaL 
edlenos 
JSESANI 
NOILVWYOSSNVYL  83IdNOs 
Wesidads GSHLOOWS NO!LINNd ONIHLOOWS ; WNYLIIdS IWNIDINO 


The truncation of the interferogram at the end 


of the sampling at re 


the interferogram by a boxcar function where the inter- 


can be thought of as multiplying 


ferogram is multiplied by one if x is less than Re 


and multiplied by zero if x is greater or equal to xe 


The shape of the truncation function is shown in Figure 
14(a). This truncation imposes the instrumental line 
shape to the spectral lines (50) of the form 


sin(21™vx,.,.)/(21vx ). The peak shape is shown in 


max 
Figure 15(a). Because of the method of computation, the 
negative side lobes of the sine x/x function appear 
positive. The peaks in Figures 15 and 17 were generated 
by Fourier transforming a 4096 point interferogram 
consisting of a cosine wave sampled at 4 points per cycle 
with the appropriate apodization. This is the resultant 
instrumental line shape when there is no intentional 
apodization, that is, when the interferogram is not 
multiplied by a function after the interferogram is 
acquired. When the term no apodization is used, it means 
that no intentional apodization has been performed on the 
interferogram. The resolution that can be obtained in 
this case is given by 1/2x_,,. For most of the spectra 
presented here, the interferograms are 4096 points long 
with the direct sampling rate and equal number of data 


points on either side of the ZPD position; therefore 


Kaax=0+1296 cm (2048 X 6.328 X 107° cm). The nominal 
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Figure 14, Apodization functions: 
(a) boxcar, (b) triangular. 
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resolution is about 4 emane However, the peak shape 
with the truncated interferogram gives large side lobes, 
the first side lobe being about 20% of the peak maximum, 
i many cases, the side lobes can cause confusion, 
especially if small peaks are close to large peaks, 

To improve the instrumental line shape, the 
interferogram can be multiplied by a mathematical function 
before the Fourier transform is performed (23,25,55-57). 
This function is usually normalized so that the value is 
Petween zero and one. In general, apodization functions 
do not necessarily have to reach zero. The apodization 
is performed easily since the interferogram is already 
stored in the computer. Typical functional dependencies 
include linear, exponential and Gaussian function. 

A consequence of apodizations that give more weight 
to small x and less to large x is the improvement in 
signal-to-noise ratio in the final spectrum, This comes 
about because the rapid changes (high spatial frequencies ) 
in the spectrum are represented in the interferogram at 
larger optical retardations than the slower changing 
information. Noise tends to be of higher spatial 
frequency than spectral peaks, therefore the noise is 
reduced, But there is the standard trade-off of 
resolution and S/N ratio. The spectral peaks also 
contain some high spatial frequencies which determine the 


sharpness of the peak, Therefore, when the noise is 
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reduced by apodization, the resolution decreases, 

The boxcar apodization function shown in Figure 14(a) 
can be implemented by a routine that multiplies the inter- 
ferogram data points less than Nl and greater than N4 by 
zero and leaving points between N2 and N3 unchanged. 

OOS linearly decreasing function, called triangular 
apodization (Figure 14(b)), can be implemented with the 
routine by setting N2 and N3 to the data point corresponding 
to the ZPD position, All points between Nl and N2, and 
between N3 and N4 are multiplied by the y value of the 
slope which varies linearly between 0 and 1. Points less 
than Nl and greater than N4 are set equal to zero. The 
instrumental line shape resulting from triangular apod- 


2 


ization is sin- (2uvx )/(2nvx 3) » The peak shape is 


max 
shown in Figure 15(b). Note that the side lobes are less 
prominent than for the boxcar apodization. The resolution 
possible with triangular apodization is 1/x,,, or half the 
resolution of a boxcar apodization. For the example given 
earlier (4096 point, double-sided interferogram), this 
means a nominal resolution of about 8 mea. 
The type of apodization found to be more useful 
with narrow emission type signals is the Gaussian function. 


The form of the function used is 


FCN(I) = exp( -FAC os Ne) Gp) 
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where I=point number, Z=point number of the ZPD position, 
NP=total number of points, and FAC=a factor that determines 
the width of the Gaussian function, The Gaussian apodization 
is shown in Figure 16. The ZPD position is shown at 
point 2048 (for 4096 point double-sided interferograms). 
Plot A results when the variable FAC=10.0, plot B when 
FAC=5.0, and plot C when FAC=2,.5, The peak shapes resulting 
from Gaussian apodization are shown in Figure 17(a) for 
FAC=2.5 and in Figure 17(b) for FAC=5.0. The resolution, 
when Gaussian apodization is applied, should be similar 
to triangular apodization if FAC is approximately 2.5. 
For larger values of FAC, the function is more heavily 
damped and therefore resolution will be lower. The side 
lobe suppression with the FAC=2.5 Gaussian apodization 
is more effective than with the triangular function while 
retaining similar resolution. The FAC=5,0 Gaussian 
function almost reaches zero at maximum optical retardation 
with the result that side lobes are not evident (actually 
they are less than 0.5% of the peak amplitude); however, 
the resolution is lower. 

A practical example of apodization is shown in 
Figures 18 and 19. The flame emission of a sodium 
solution was measured (details in Chapter VII). Figure 18 
shows 3000 points of the 4096 point double-sided inter- 
ferogram that resulted from time averaging 50 successive 


scans, That two emission lines are present can be deduced 
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Aor 4096 
POINT NUMBER 


Gaussian apodization functions: 
A--FAC=10.0, B--FAC=5.0, C=--FAC=2,5. 
See text for details. 
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Peak shapes resulting from Gaussian 
_apodization: (a) TFACS2,5,° (b)PrAG=5. 0. 
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Interferogram of Na 589 nm doublet (3000 of the 4096 points shown). 
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Interferogram of Na589 nm doublet after Gaussian apodization (FAC=2.5). 
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from the interferogram. The beat pattern is the result of 
the summation of two cosine waves of slightly differing 
frequencies, Figure 19 shows the result of multiplying 
the interferogram shown in Figure 18 with a Gaussian 
function (FAC=2,5). Again only the middle 3000 points 
of the 4096 point interferogram are plotted. The 
Fourier transform is performed on both interferograms and 
the amplitude spectra are shown in Figures 20(a) and 20(b). 
Only the part of the spectrum around the Na doublet is 
plotted to show the instrumental line shape that has been 
imposed by the apodization. It can be seen that the 
spectrum from the apodized interferogram is more amenable 
to interpretation, 

The apodization process does not give any more 
information than is present in the unapodized spectrum; 
in fact, information can be lost by apodization. This 
can occur if the apodization results in a lower resolution 
and causes two lines to become unresolved. Therefore 
the implementation of apodization must be done with care 
to give the desired instrumental line shape without 


removing essential information in the spectrum. 


60 


vo een & wits ze: nie tae 


whos tee vitoall? 


« : 
<> 
4 


7. bbe Yas fipeaayt 


ack hes4giq SIH fer eats ine 


‘te eth 


cae 
‘ye ae yleey shah r 


tJ mien tts 2X8, since 
brits tye ae ae ua a 
mit Peddemessaes* Ad wecuty on 


ano dt ete tboge ua. 
s - 
ik, blot biegtey amy satel 
i f os. 
an | 
j © j hh (a 3 ; £4 
| he 


UL Lu 
— OQ 
—) a, 
face b~ 
a a 
=. = 
a < 
(a) (b) 
omen a es 
10880 17080 16880 17080 
ee -1 
cm cm 
Figure 20, Na 589 nm doublet spectral peak shape: 


(a) with no apodization, (b) with Gaussian 
apodization (FAC=2.5). 
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CHAPTER V 
PHASE PROBLEMS IN FOURIER TRANSFORM SPECTROSCOPY 
—————— aR EU DEEGLERUOGOPY 


The technique of Fourier transform spectroscopy 
involves the recording of the interferogram that results 
from the two beam interference of light via a Michelson 
interferometer. Ideally, the interferogram is a 
symmetric function of the optical path difference. 

That is, the function is the same from 0 to Teas 


ax 
from 0 to =x - Therefore, only one side need be 


max 
recorded in this ideal case, 

However, several factors contribute to give an 
asymmetric interferogram, If all the optical and elec- 
tronic components were perfect, at zero path difference 
(ZPD) every spectral frequency would be at zero phase, 

But this is not the case and phase errors are always 
present, The fact that the light passes through real and 
not perfect optics may cause different spectral frequencies 
to have slightly differing phases at the ZPD position. 

Also the sampled interferogram may not contain a data 
point exactly at the ZPD position. This results in a 
constant phase error due to sampling. Detector electronics 
may have phase response that is dependent on the frequency 
of the input. This is especially true if electronic 


filters are used to limit the frequency bandwidth of the 


signal. Hence, imperfect optics and measurement 


62 


a age he 
; ale 

ha sari: aq' ee >t aries aa va Shir as eee ees iF 
1 SBA th. re na “a5: af ang’ Ve giibrobas 8 aN a 
ie toes ec ele opal'T 3a scheint capt get 


a Se che yonGd eam! a atts ci ai a FIM 7 


Brystky Wray tap eee oft Ww nity acl 
‘ae: alee vata at. @ cosas oor 


eof SA. et 7 Lito, yeti wer? ‘anal , 
ok 4 Saetel ele oe 
= my 
10 | 4 eee 5 satan eee 


® 


sity ike? 4a Desitadtetaeted 
ce a ey shape bi, *rithta’ Xi 
ee fai hia tana oda 


; AiG Dart Shae ete ro ee 


A ' f if 
Fi ies ita. $278 ee Ag 2 sat a 


r 1%) Aha e ‘ y sta i ih XG a ebanieit i v 


4 ' 
rg 


; shia } pL ei ou 
/ a7 ie feta ag A ia ngls 38 qa 
PU a ubidves es by Tera Det: 
pen) GF 0, Jehovah AC 44 Senger bap 


ws § , Vid Tess a bay 


oe 


Asie. ew a ¥ snare 4 a Lents Jd are 
ae! 


“ig 


earn i pei BAe §1Ig) (oars sents 


63 


electronics produce phase errors that are varying functions 
of the spectral frequencies, 

The detailed mathematical description of phase 
errors and correction methods have been well documented 
since it is a problem that is always present with single- 
sided interferograms and must be corrected to give an ac- 
curate output spectrum (15,23,58-61). The discussion 
given here will be of a qualitative nature and will 
emphasize the phase problem with respect to line emission 
type signals in the UV-visible-near IR regions. 

A special problem arises when the interferogram is 
sampled at a rate that results in aliasing. Since the 
phase error is dependent on the spectral frequency of the 
Signal, each aliased region requires a different phase 
correction, This makes it next to impossible to phase 
correct an aliased spectrum. It should be mentioned 
that problems with phase correction of line emission 
spectra in the IR ( not aliased) have been reported 
(62563). 

Measurements taken with the first data system 
using the PDP-8/e minicomputer will illustrate the 
phase problem. Only single-sided interferograms could 
be acquired at that time, The flame emission Ot eDike Ks 
Rb, and Cs was studied. Figure 21 shows the spectrum 
from a potassium solution (1000 ppm). Spectrum A 


results from the Fourier transform of the triangularly 
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Figure 21, Potassium spectrum from a single-sided 
interferogram (512 points). 
A: no phase correction, B: phase 
corrected. 
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apodized interferogram. The K doublet is obviously 
very distorted, 

Phase correction procedures utilized in Fourier 
transform NMR (64) were attempted on these flame emission 
spectra. The procedure for phase correction involves 
convolving all the points in the real and imaginary 
Spectra with a cosine and sine wave. The phase angle 
required to correct a peak is obtained by inspection of 
the line shape in the real or imaginary spectra. If the 
spectral peak needs no phase correction, the real and 
imaginary outputs from the FFT would resemble Figure 22(a). 
If the real part were distorted, the phase angle is 
estimated from an inspection ef a chart of peaks distorted 
from O through 360 degrees, Figure 22 shows peaks at 
0,30,90,120, and 180 degrees of phase error. 

The phase corrected K spectrum is shown in spectrum 
B of Figure 21, Considerable improvement has been made. 
Similar phase correction procedures were performed for 
the Li, Rb, and Cs spectra and the results are shown 
in Figure 23, The analog interferograms of the four 
elements are shown in Figure 24. The digitized inter- 
ferograms sampled at 2.5312 um intervals are shown in 
Figure 25, The aliasing can be seen in the digitized 
interferograms, If sampling was done at a rate high 
enough to avoid aliasing, the digital interferogram 


would have the same appearance as the analog signals. 
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Figure 22, 
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Real and imaginary output peak shapes for 
phase errors of 0, 30,90,120, and 180°. 


(a) 


(b) 


Cc) 


(d) 


(e) 


66 


- 
-_ 
_ 
— 


sg dv ee haitien, Arh ee oeoA 
7 a hd ae Ni j a) % . Mi 4 ie | 0 : a.) Ca 15 o 


67 


°sp pue ‘qy Sy ‘TT Jo earqoeds peqoezi00 eseug °¢z7 eANSTy 


i- 


429796 hynileaes ISetl 4? 9z8"E1 Be ay Bata Amey gle ae 


V £768 8D | 
| 
ly cre a 
vizseso| G vooez qa; | 2) 
1? 9eBEL wnissv1Od sell, M>-978'E! WNIHLN Coss 
ee ee ee ee en es Big ens elle ee x ane en) 
ee ec ee Cae ee 
| 
| | 
: 
i 
| i 
v 8692 > 
| 
V 799» | 9 v 202911 | y 


POTASSIUM 


LITHIUM 


CESIUM 


RUBIDIUM 


Analog interferograms of Li, K, Rb, and Cs. 


Figure 24, 
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However, when a solution of all four elements was 
measured (Figure 26), the spectrum was very distorted and 
similar phase correction procedures could not improve the 
spectrum to any great extent. This is because the spectral 
lines are in different aliased regions and each region 
requires a different phase correction angle. The cumulative 
effects of the use of different correction angles for each 
spectral region produces distortions of its own,thus 
making overall correction of the spectrum impossible. 

A method of phase correction that has been fairly 
successful in Fourier transform infrared measurements does 
not seem to work as well for line emission spectra. The 
reason it works for most IR measurements is due to the 
fact that these spectra are of the absorption type; that 
is, there is broadband information present. This allows 
the calculation of a reference phase spectrum which can 
be used to correct the sample spectrum. A low resolution 
reference phase spectrum is satisfactory and this can be 
Calculated from a small portion of the interferogram on 
both sides of the ZPD position. The phase correction 
procedure is described in reference (23). For line 
emission type signals, there is little or no broadband 
information; hence, an accurate reference phase spectrum 


is difficult to obtain. This makes phase 
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Figure 26, Multielement spectrum of Li, K, Rb, and 
Cs (not phase corrected). 
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correction of line emission type signals not very 
successful, 

The problem of phase correction can be totally 
eliminated if the interferogram is sampled from -x 


max 


to +x 


max? that is, a double-sided interferogram is 


acquired. If the amplitude spectrum is calculated 
(the amplitude spectrum is the square root of the sum 
of the real and imaginary parts squared), the results 
are independent of the phase (Chapter 12 of reference 32), 

This is shown in Figure 27, A double-sided 
interferogram was acquired using a hollow cathode lamp 
as a source (most of the lines shown are due to the 
neon gas filler of the lamp). Triangular apodization 
was applied and the result is spectrum (b). The peaks 
show no distortion. To produce the single-sided inter- 
ferogram, the points to the left of the ZPD position 
were set equal to zero. The apodizing functions are 
shown above the spectra. Triangular apodization was 
applied and the result is shown in spectrum (a). 
Distortions similar to the flame emission data (Figure 26) 
can be seen, 

The flame emission measurements of Li, K, Rb, and 
Cs were taken using double-sided interferograms and the 
results are shown in Figure 12(d) in Chapter III. This 
spectrum and the one in Figure 26 can be compared 


directly. It should be noted that the improvement in 
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quality of the spectrum in Figure 12(d) occurs in spite 
of the fact that resolution is lower by a factor of two 
(same number of data points but a single-sided inter- 
ferogram samples to twice the optical retardation of a 
double-sided interferogram), 

| The question arises as to why double-sided inter- 
ferograms are not used exclusively, There are several 
reasons. A major reason is that double-sided interferograms 
need twice as many data points for the same resolution 
compared to single-sided interferograms, The resolution 
depends on the maximum optical retardation, x 


max ° 


Sampling from -x 


+ i he same resolution 
Fae to Ke auetves te 


as from 0 to va With twice as many points the 
computing time for the Fourier transformation is greatly 
increased, Also the computer memory needed must be 
doubled if an "in-core'' FFT is used, Other factors are 
that the time required to acquire the double-sided data 
is doubled, and the moving mirror must be driven twice 
as far, 

The above considerations make the acquisition of 
double-sided interferograms impractical for extremely 
high resolution measurements where in the order of 
hundreds of thousands of points need to be acquired. 

However, for line emission type signals and medium 
resolution, the double-sided interferogram is not 
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are to be obtained, because of the lack of adequate phase 


correction procedures for line emission signals. 
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CHAPTER VI 
THE MICHELSON INTERFEROMETER SYSTEM 


A, Description of the Michelson Interferometer 


The current Michelson interferometer has evolved 
over several years, The initial interferometer design 
was developed at the Department of Aerospace Engineering, 
High Altitude Engineering Laboratory at the University 
of Michigan by L. W. Chaney in cooperation with Goddard 
Space Flight Center for balloon borne and space 
experiments (65,66). A non-repetitive slow scanning 
mirror drive system was built at the University of 
Illinois (52). It has been substantially modified and 
developed into the current interferometer design. This 
chapter will describe the interferometer in detail. 

A block diagram of the interferometer system is 
shown in Figure 28, There are three optical inputs to 
the interferometer: a He-Ne laser, a white light source 
(tungsten bulb), and the spectral signal of interest. 
Laser fringe referencing is used to sequence digitization. 
The laser fringes are also used to control and stabilize 
the velocity of the moving mirror using phase locked 
loop techniques. The mirror drive system consists of an 
electromechanically driven piston supported by an air 


bearing and a mirror mounted on the piston. The mirror 
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Figure 28, Block diagram of the Michelson 
interferometer system. 
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movement is repetitive and both the scan rate and length 
can be easily set or altered, With the control signals 
derived from the white light interferogram and the laser 
fringes, spectral signal interferograms can be precisely 
time averaged, A unique "pretrigger" approach allows 
acquisition of any desired number of data points before 
and after the zero path difference position. All three 
optical inputs share the same beamsplitter and Michelson 
mirrors. This feature greatly facilitates alignment of 
the interferometer, 

The concept of the design is to be flexible in 
order to ease alignment, assembly, and modifications. 
Therefore the interferometer consists of modules, each 
of which can be changed easily. For example, the source, 
beamsplitter, and detector can be changed and the inter- 
ferometer realigned in a few minutes, A schematic 
diagram of the interferometer is shown in Figure 29, 

It will be useful to refer to Figures 28-30 and the 
photographs shown in Figures 31 and 32 during the 


Following discussion. 


1) Central Mounting Cube 


The mounting cube is the central unit of 
the interferometer. All the other modules mount onto 
the cube. The cube is approximately 10.8 cm (4.25 inches) 
and was machined out of a solid piece of aluminum. The 


entrance port visible in Figure 31(a) is 6.05 cm (2.38; ing) 
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DETECTOR 
MOUNT 


(b) WHITE LIGHT 
REFERENCE ASSEMBLY 


Figure 30, ng Lens detector assembly. 
b) White light reference assembly 
schematic. 7 
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Figure 31, 


Photographs of the interferometer: 
(a) front view, (b) top view. 
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Ficure 32, 


Photographs of the interferometer: 
a) mounting cube partially disassembled, 
b) beamsplitter and white light assembly. 
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Fieure 32 (c). 


The interferometer 
disassembled, 
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in diameter. In order to allow the He-Ne laser input 
module to be mounted from the bottom, the cube is mounted 
on a brass plate over a hole in the optical table, 

2) Beamsplitter Assembly 

The beamsplitter assembly is best seen in 

Figure 32(b) where it is raised out of the central cube, 
It is mounted under the top face plate of the interferometer, 
The assembly consists of two optical flats (a beamsplitter 
substrate upon which the beamsplitter coating is deposited 
and a compensator) mounted on two aluminum holders 
separated by an aluminum spacer. The beamsplitter 
substrate and compensator are each 7.62 cm (3.00 in.) in 
diameter and 0.953 cm (0.375 in.) thick and their dimen- 
sions should be matched as closely as possible. For 
work in the UV, visible and near-IR, fused silica flats 
were used (from L. H. Sampson Co., 22730 Orchard Lake 
Road, Farmington, Michigan 48024) with flatness to %/10 
in the visible, For work in the mid-IR, NaCl flats were 
used, The salt flats are flat to better than 1A in the 
visible (obtained from Frank Cooke Inc., 59 Summer Street, 
North Brookfield, Mass. 01535). 

Germanium was used as the beamsplitter material 
and was coated on the substrates by a local optical 
shop at the University of Alberta. The germanium was 
coated on until the transmittance was approximately 50% 


in the spectral region of interest. Figure 33 shows the 
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Figure 33. Transmission spectra of the beamsplitters: 


a) quartz substrate for UV,visible, 
b) NaCl substrate for mid-IR. 
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transmission spectra (light incident at 90°, air 
reference) of the fused silica substrate and the NaCl 
substrate beamsplitter assemblies, 
3) Fixed Mirror Assembly 

The fixed mirror assembly consists of a 
spring loaded ball-in-socket mount with two sets of 
opposing spring loaded plungers and screw-driven 
adjustments. This design has proven to be very stable. 
The fixed and moving mirrors are both 4.2 cm (1.64 in.) 
in diameter, 0.64 cm (0.25 in.) thick, flat to A/10 
in the visible and coated with aluminum. (The mirrors 
were obtained from L. H. Sampson Co.). 

4) Moving Mirror Assembly 
The moving mirror assembly is based on an 

air bearing suspension system driven by a standard loud- 
speaker magnet and coil assembly. (The air bearing, 
Model L.B.-10, was obtained from Bearing Specialty Co., 
90 Windom Street, Boston, Mass. 02134). The moving 
mirror is mounted on a piston inserted inside the air 
bearing. The loudspeaker coil is mounted at the other 
end of the piston and rides into the magnet assembly. 
Weak retaining springs keep the piston assembly roughly 
centered in the air bearing. Tank nitrogen (20-25 psi) 
is used to float the moving mirror piston. The complete 
air bearing assembly can be shifted back and forth in 


the moving mirror assembly for coarse positioning of the 
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moving mirror static position, The magnet holder is the 
only unit of the modules not physically attached to the 
interferometer cube. It is mounted on the optical table, 
A plexiglass dust cover encloses the gap between the magnet 
holder and the moving mirror assembly. With the present 
design the maximum mirror movement is about 1 cil. | 
5) White Light Reference System 

The white light reference system consists 
of two assemblies mounted on the top face plate of the 
interferometer, One assembly contains a tiny tungsten 
bulb and collimating lens (see Figure 30(b)). The 
collimated white light source is reflected into the inter- 
ferometer by a small diagonal mirror, such that only a 
small part of the top of the Michelson mirror is used 
for this signal, At the output end, the white light 
signal is reflected back up to a silicon photodiode 
detector by another snall diagonal mirror and focusing 
lens, The arrangement of this important periscope 
system is more clearly shown in Figure 32(b) where the 
pee face plate, with attached beamsplitter has been pulled 
out of the central mounting cube, This photograph also 
emphasizes the modular nature of the interferometer and 
the fact that the beamsplitter can be exchanged easily. 

6) He-Ne Laser Reference System 
The He-Ne laser reference system is very 


similar to the white light reference system and is mounted 
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Ga the bottom face plate of the interferometer. The two 
periscopes project through the hole in the optical table 
over which the central cube is mounted, The beam from 
the He-Ne laser (Spectra Physics Model 133), mounted 
under the top surface of the optical table, is reflected 
up into the periscope and through the interferometer by 
two small right angle prisms, The laser path through 
the interferometer uses only a small part of the bottom 
of the Michelson mirrors. Since the Laser beam is 
inherently collimated, no lenses are needed in either 
periscope of the laser system, At the output periscope 
the laser signal is reflected down a custom made inter- 
ference filter to a Si photodiode detector, (The filter 
bandpass is centered at 632.8 nm and was supplied by 
imtrared Industries, Inc.,. Thin Film Products, 84 Fourth 
Ave., Waltham, Mass, 02154), The filter is mounted in 
the same manner as the lenses for the white light assembly. 
The laser periscopes can be seen in Figures 31(b)and 32(a). 
7) Detector Assembly 

Two detector assemblies have been used 
with the interferometer. One is based on a lens to focus 
the light to the detector (Figure 30(a). Both glass and 
salt lenses have been used. The salt lens was ground 
in our own laboratory from a NaCl cylinder. Both lenses 
Peveeaudiameter and focal length 0f.5.97 cm:(2s35 ins). 


However, at the short wavelengths of the UV-visible 
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region, a short focal length single lens is prone to 
Peeees tions that could affect the quality of the output 
signal. Therefore, the second detector assembly consists 
of an off-axis parabolic mirror (Figure 29, 31(b)), 
folding mirror and detector mount. The off-axis para- 
bolic mirror is mounted on a sub-assembly identical to 
the fixed mirror assembly, The focal length of the off- 
axis parabolic mirror is 16.8 cm (6.6 in.) with the focal 
point 6.1 cm (2.4 in.) from the inner edge. (The para- 
bolic mirror was obtained from Optronics International 
Sales Corp., 7 Stuart Rd., Chelmsford, Mass. 01824), 

Several detectors have been used: TGS (triglycine 
sulphate) pyroelectric for mid-IR measurements (series 
T-300 unit from Barnes Engineering Co., 30 Commerce Rd., 
Stamford, Conn. 06902), Si photocell for near-IR and 
visible, 1P28 (RCA) photomultiplier tube (PMT) for 
visible, and R166 solar blind (Hamamatsu) PMT for 
ultraviolet measurements. 

It should be pointed out that the focal length of 
the off-axis parabolic mirror represents a compromise. 
The radius of the first fringe y in the focal plane of 
the output of a Michelson interferometer is equal to 
£(A/a) 1/2 where £ is the focal length of the output mirror, 
A the wavelength and d the optical retardation. To get 


maximum sensitivity, it is desirable for the detector to 


see the entire central (zero) fringe. A smaller detector 
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aperture will give a lower signal level. But an aperture 
larger than the central fringe will degrade the signal 
because the detector will be observing more than one 
fringe simultaneously. 

Thus, for optimal IR performance in the 2-L5um 
région, the focal length of 16.8 cm.is too long as it 
results in a fringe size (y=1.47 cm for 10um wavelength 
and d=0,13 cm) that overfills our TGS detector (element 


1 mm x Ll mm) by a considerable amount, particularly at 


the longer wavelengths, This results in lower sensitivity. 


However, for visible and especially UV performance, the 
longer focal length is desirable. Typical IR detectors 
have small sensor areas and thus act as their own aper- 
tures. With a PMT, a circular aperture must be mounted 
in the exit focal plane in front of the PMT. A short 
focal length, combined with working at wavelengths over 
an order of magnitude shorter, would result in a very 
small fringe radius making fabrication of an appropriately 
small aperture difficult and alse making focusing much 
more critical, For example, with f=5.08 cm (2.0 in.), 
250 nm wavelength and d=0.13 cm, y=0.07 cm. Also, the 
problem becomes increasingly severe as one goes to 
higher resolution measurements; that is, as d increases, 
y decreases further. Since our initial experimental 
goals were to investigate the applicability of our 


Fourier transform spectrometer to UV and visible 
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measurements, we chose a longer focal length (16.8 cm) 
than would be optimal for a purely IR instrument in 
order to use larger exit apertures with the PMT detectors, 
Because of the modular nature of the interferometer, it 
would be easy to design a separate detector assembly 
optimized for IR performance, 
8) Input Optics 

The spectral input to the interferometer 
is collimated by an off-axis parabolic mirror (identical 
to the one in the detector assembly). The source is 
placed at the focal point of the mirror. Sources that 
have been used are hollow cathode lamps and a modified 
Varian burner (see Chapter VIIA for details) for UV- 
visible-near IR measurements and a rapes for mid-IR 
measurements. A filter holder for 2 in. x 2 in. filters 
can be placed at the interferometer entrance port if it 
is necessary to block unwanted portions of the spectral 


region. 
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B, Consequence of Using Only One Interferometer 


The design of our interferometer uses the same 
Michelson mirrors and beamsplitter for all three optical 
signals, The white light channel uses the top section 
of the Michelson mirrors, the He-Ne laser channel the 
bottom section and the spectral signal the center of the 
mirrors. This is a major difference compared to most 
commercially available systems which use two or even 
three independent interferometers for these signals. 

We have found the single interferometer design to be 

both simple and effective, In particular, alignment and 
the maintenance of alignment are greatly facilitated since 
only one mirror needs to be adjusted. 

If the interferometer is grossly misaligned, as it 
would be if the beamsplitter were changed, the procedure 
for alignment is as follows. The detector assembly is 
easily removed and, using the laser signal, coarse 
adjustment is made visually by lining up the laser beams 
from the two arms of the interferometer. Then, with the 
moving mirror scanning, the fixed mirror is fine tuned 
to locate and maximize the laser fringes signal. The 
alignment at this point is nearly optimal. The white 
light interferogram signal is then maximized. With these 
two signals aligned, the spectral signal channel is 


automatically aligned as only one mirror is involved. 


week ats a Sesotho 
har’ ‘ ae ee 
| Hh 
i Ly 1. : re 
‘ihe Sr RRS Sa Seni ee ue, 
Li : ae valid ; 
faetsqe see Pia Sin SEE | 


ane Sa rea Lem 


* s pre hth led 
Si 2 Shrreio oSae F cay yt 
' ’ aif : 5 ya i j 7 ; « fe i . 
3h)~ tf Cle = ff ~ I hss C ‘% hi ? oy : 5 . ‘¥ ; . sore sas y ay 
Be , Pi Be ie i. a) 5989) {os ne bs cain a, aa, ban 


yoys mW OW | ey et, ‘ it? "ES, We ren q ue ‘seat oni 


 » _ . | i 7 


JT aaot | pee dione 80% ca cell dani og ae . 


Lee te nds re, 


has Clee 2 See iva 


ibe AO meh rion ke (ate at | 
we 


a 


Pitas? WIIee te Sse) eee ye 20 sane 
Say Sy as | A WAAR ROR ‘shor 


Yair srs Wag td 


» 
w 
Le 


aH hal Sh ee pad ome Bs Ae hes pieeeaKy eit! my seal 


: 
j 


ai 4 


geen) 3 Sis ‘ay 1 i eld, orn a ised ‘ ra pevorme 4 

{ _ 

a9 pene? 

2 ne Ets | 1ae8 bf ei (Pe) Lt oh rhige cis "1 meee ay ‘mae at 7 a 
babs : ag 

and, matw _ster CRUAL is ok kell Shree y vet 2 + RPS: ane “eo i 


< 


ve 
' 
iT 


hans  octa eF Texan Hart? aia getisagnad ‘eras a 


Je 
Af 


at, beaels Gegubti “sont ets Suntaan aa ey 
wate eet  laeicon elspen el aatag etrks a0, 7 


{ 
, 


- 
ry 
Lid 
a 
_ 

i] 


mevis Ash. ~boslueieen ried et) DRaR ee ae goux08 


af fagnetA tahgte Inadeae wee vera t Tp ats 


we 


cbeviewnt wh cessin she. TIA (a Sha Lie: 


= 


es) 


In contrast, systems with multiple interferometers can 
be very time consuming to adjust; particularly the spectral 
Signal and white light channels, 

Since the thickness of Ge required for an effective 
mid-IR beamsplitter is such as to be essentially opaque 
in the visible, this design does require that small 
windows of higher transparency be left on the top and 
bottom of the salt beamsplitter to allow the white light 
and laser signals to pass, This is easily achieved via 
a two step coating procedure using a mask on the second 
coating. In addition, even though salt flats can only be 
obtained with a one wavelength flatness specification, 
the white light and laser fringes are of excellent quality 
with our mid-IR salt beamsplitter assembly. For the 
near-IR and visible beamsplitter, one coating is 
sufficient; windows of different “ZT being unnecessary 
for the white light and laser. 

Another major consequence of the single inter- 
ferometer design is the coincidence in time (and mirror 
position) of the white light and spectral signal inter- 
ferograms, The occurence of the white light interferogram 
peak is normally used as a precise start flag when time 
averaging repetitive interferograms, This means that 
only single-sided interferograms can be acquired with 
no points acquired before the central fringe Oletne 


spectral signal interferogram. Single-sided interferograms 
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make it difficult, if not impossible, to phase correct 
the spectra properly as shown in the chapter dealing 
with phase. In fact, this is the major reason why, in 
most designs, a separate interferometer is used for the 
white light channel. Thus, one of the Michelson mirrors 
of the white light channel can be offset so that the 
white light central fringe occurs several hundred laser 
fringes before the spectral signal central fringe. That 


is, the zero path difference(ZPD) point for the white 


light channel occurs before the ZPD point for the spectral 


signal, 

To overcome the single-sided interferogram problem, 
a versatile hardware and software system was developed to 
precisely acquire and time average any desired number of 
data points both before and after the signal central 
fringe. The procedure involves generating a start flag 
(some imprecise point before the ZPD position) to 
initiate taking data. Then a stop flag (precise because 
it is generated by a preset number of laser fringes after 
the white light central fringe) is used to halt the data 
acquisition. The time averaging is then performed from 
the last data point backwards. This "pretrigger" 
approach and how it is accomplished is explained in 


detail in the next section. 
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C. Interferometer Electronics 


Uhe interferometer control electronics is shown 
in Figure 34 and the timing sequence of the various 
electronic signals is illustrated schematically in Figure 
35. The two key control signals generated by the inter- 
ferometer are shown in Figure 36(a). The white light 
fringes signal measured at the output of OAl is shown 
in Figure 36(a) (upper trace). The laser fringes signal, 
simultaneously measured at the output of OA2, is shown 
in Figure 36(a) (lower trace). These signals will be 
referred to as the white light analog and laser analog 
Signals, The white light signal, in combination with the 
laser signal, is used to generate a precise stop flag in 
order to sequence averaging of many scans. The laser 
eienal is used to stabilize the velocity of the moving 
mirror and to initiate digitization accurately. 

1) Moving Mirror Drive 

The basic mirror drive signal is a tri- 
angular waveform obtained from a function generator 
(Heath Model EU-81A). This signal is amplified by a 
power amplifier (Hewlett-Packard Harrison 6824A Power 
Supply Amplifier) to provide the current to drive the 
loudspeaker coil. The drive rate and length are set 
with these two units. Most of our measurements have 


been carried out with a drive frequency of about 0.2 Hz 
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Figure 


Sie 


(a) 


(a) White light analog (upper trace)and 
laser analog (lower trace) signals. 

(b) White light analog (upper) and white 
light digital (lower) signals. 

(Scope time base: 0.2msec/division) 
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Bieure (37, 


(b) 


(a) Phase locked loop control of laser analog 
Signal: upper, no PLL control; lower, 
WLEN PLE coneEro.. 

(b) White light analog (upper) and laser 


digital (lower) signals. 
(Scope time base: 0.05 msec/division) 
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and an amplitude sufficient to drive the mirror 0.3-0.5 
em, This produces laser fringes at a frequency around 
4kHz. 

2) Laser Channel 

The signal interferogram must be sampled 
at accurate distance intervals, For each cycle of the 
laser fringe signal shown in Figure 36(a) (lower trace) 
the mirror moves 0.3164 um and the optical path length 
changes 0.6328 um. Thus this signal can be used to provide 
an accurate clock for digitization. The frequency of 
this clock is determined by the velocity of the moving 
mirror, Even though this clock will accurately mimic 
any velocity changes in the mirror, it is still desirable 
to stabilize the mirror velocity because of the variable 
effects that low and high pass electronic filters will 
have on an interferogram signal whose frequency composi- 
tion is changing due to mirror velocity changes. In 
our system, the mirror velocity is stabilized using a 
phase locked loop feedback control system. 

A phase locked loop (PLL) consists of a phase 
comparator, a low pass filter and a voltage controlled 
oscillator (VCO) (54,67) and it can be obtained as a 
single integrated circuit (IC) package. In our circuit 
we have used a Signetics 565 PLL, Functionally the PLL 
compares the phase and frequency of an imput frequency 


to the internal VCO frequency and generates an error 
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voltage proportional in sign and magnitude to the phase 
and frequency difference of the input and VCO frequencies, 
This error signal is applied through the low pass filter 
to the VCO control input (internal connection on the iG) 
and drives the VCO until it is at exactly the same 
frequency as the input frequency, that is, the VCO output 
becomes locked to the input frequency. Clearly, no 
stabilization of the mirror velocity is achieved merely 
by just locking the VCO of the PLL to the laser frequency. 
But the error signal discussed above is available as an 
output from the integrated circuit (pin 7). This error 
signal is applied, with appropriate polarity, to the 
external voltage control input of the mirror drive 
function generator. 

This input of the function generator allows an 
offset to the basic drive frequency (about 0.2 Hz) so 
that the ae euat drive frequency depends on the polarity 
and amplitude of the error signal. For example, if the 
error signal is momentarily positive, the drive frequency 
is greater than 0.2 Hz for that time. This feedback 
loop (OA3 and OA4) stabilizes the mirror velocity and 
hence the laser fringe frequency. The capability of 
this system in stabilizing the mirror velocity is shown 
in Figure 37(a). The upper trace is the laser analog 
signal without feedback control and the lower trace 


with feedback control. Each trace is a 0.5 second time 
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exposure taken during one scan of the mirror drive, 
Complete details on the operation of the 565 PLL 
can be obtained from Signetics (68). Only those 
connections pertinent to this application will be discussed. 
The VCO free running frequency is set by the value of the 
capacitor to ground from pin 9 and the resistance to +1L5V 
from pin 8. The frequency is set by the 10k ohm 
potentiometer to about 4 kHz. The 0.5 uf capacitor from 
pin 7 to +15V sets the internal loop low pass filtering 
and the 0.001 uf capacitor between pins 7 and 8 is for 
decoupling. Pin 7 is the error voltage output and herice 
provides the feedback signal for the mirror drive 
stabilization loop. Pins 2 and 3 are the external inputs 
to the phase comparator. The laser analog signal is 
connected to pin 2. The voltage divider network on pins 
2 and 3 is necessary for peoner operation of the PLL 
integrated circuit from a single power supply. Pin 4 
is the VCO output which, when the complete system is 
interconnected, is phase locked to the stabilized laser 
frequency. Pin 4 must also be connected to pin S,awhicn 
is the internal input to the phase comparator in order 
to close the loop on the integrated circuit. To complete 
the laser channel circuit, the VCO output is converted to 
a, iTL Tevet signal by the OA comparator. The output of 
the comparator will be referred to as the laser digital 


signal (Figure 37(b), lower trace) and is used to clock 
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the analog-to-digital converter (ADC). 
3) White Light Channel 

In order to time average interferograms, 
some type of precise control signal is required to insure 
coherent addition of repetitive scans. This is the primary 
function of the white light interferogram. As discussed 
in Chapter VIB, in the single interferometer design, the 
white light and signal interferogram central fringe occurs 
coincidentally, Since the white light fringe is the only 
accurately known point (zero path difference position) 
during the mirror scan, using it to initiate data 
acquisition results in single-sided interferograms (no 
data points acquired before the ZPD position of the 
signal interferogram). Chapter V discussed the problems 
with calculating acceptable spectra from single-sided 
interferograms, 

Our approach to taking data 4: both sides of the 

ZPD position is to use the white light fringe to generate 
a precise stop flag. A flag is generated at some 
imprecise point before the ZPD position to start acquiring 
data points. When the white light fringe occurs, it is 
used to precisely start a modulo=N counter that counts 
the laser digital clock, The counter overflow pulse gen- 
erates a stop flag which is used to terminate data 
acquisition. Repetitive interferograms can be coherently 


time averaged by sequencing addition from the last point 
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acquired rather than the first, in a sense, time averaging 
backwards, This is achieved easily with software. This 
system will now be described in detail, It will be 
helpful to refer to Figures 34 and 35 in this discussion. 
The imprecise start flag is generated from a 
"sync'' signal available from the function generator. It 
is a square wave of the same frequency as that of the 
mirror drive triangular wave. The positive edge of the 
sync. signal, which is inverted by NAND-2, triggers 
MONO-1 (monostable) which in turn triggers flip-flop-1 
(FF-1) generating the imprecise start flag (the negative 
going edge of NAND-4),. This initiates data acquisition. 
The approximate position of the start flag, with respect 
to the white light fringes, can be set easily using 
MONO=1, It is normally set to occur about 2200-2400 
laser clock pulses before the ZPD position. The laser 
clock pulses are counted with a Heath Model EU-805 
Universal Digital Instrument in the gated counter mode. 
The clock signal is applied to the counter input and the 
counter gate is a pulse derived from the exclusive-or 
(XOR) of NAND-3 and NAND4 signals. Hence the counter 
registers the number of clock pulses between the start 
and stop flags. The modulo-N is 2048 so the total count 
is adjusted to read approximately between 4200 and 4500, 
that is, 100 to 400 more points than necessary to allow 


for drift and jitter. The start flag is imprecise with 
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respect to the ZPD position because it is essentially 
impossible to oscillate an electromechanical assembly such 
as the moving mirror assembly and have it turn around at 
interferometrically accurate amplitude. The jitter in the 
start flag has a standard deviation of about 9 laser clock 
fringes, 

The syne. signal is used for one other key function. 
Data is acquired on only one half cycle of the moving 
mirror oscillation. White light fringes are generated on 
both half cycles. The sync. signal is connected to NAND-1 
to inhibit the generation of a spurious stop flag on the 
other half cycle of the moving mirror oscillation. 

When the white light interferogram occurs, it is 
amplified by OA-1 then converted to a TTL level digital 
signal by a comparator, To provide noise immunity, the 
trigger level is gee so that only the central three or 
four fringes fire the comparator. The white light analog 
and digital signals are shown in Figure 36(b) upper and 
lower trace respectively. The first pulse of the white 
light digital signal triggers MONO-2 which immediately 
triggers FF-2, the modulo-N counter gate. This allows 
the counter (N=2048 for our present interferogram length 
of 4096 points) to count the laser digital clock pulses. 
MONO-2 pulse width is set wide enough to block out sub- 
sequent pulses of the white light digital signal. The 


overflow (carry) pulse of the counter triggers FF-3 
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meenecatinge the stop flag which is precise in terms of 
mirror movement with respect to the white light fringe ; 
(ZPD position). This stop flag (negative going edge of 
NAND-4) terminates data acquisition by the computer. 
Time averaging of the acquired interferograms in the 
computer is sequenced backwards from the last data point, 
thereby providing coherent addition of interferograms. 

Flip-flop 3, when it is triggered, also clears FF-2 
thus closing the counter gate and preventing multiple stop 
flags. The second half cycle of sync. pulse clears FF-1 
and FF=3 to be ready for the next scan sequence. The 
purpose of NAND=-3 and NAND=-4 is to buffer the flip-flops 
to meet the input characteristics of the data acquisition 
system, 

The first pulse of the white light digital signal 
must be reproducible from scan to scan. These signals 
have been observed over several hundred scans. Also, over 
500 repetitive interferograms have been time averaged and 
there has been no evidence of jitter in the first pulse 
of the white light digital signal. 

4) Detector Electronics 

The detector used for IR measurements was a 
TGS (triglycine sulphate) pyroelectric detector. A CrrCunt 
was provided by the supplier to interface to the detector 
to provide a voltage of low output impedance. This voltage 


output was connected to a low noise preamplifier (Princeton 
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Applied Research Model 225), The input is ac coupled and 
the preamp provides variable low and high pass filtering 
and variable gain. The bandpass was set at 300 Hz to 10 
kHz. The preamp output is the signal that is applied to 
the analog-to-digital converter (ADC). 

A Si photodiode was used for visible and near-IR 
measurements, The output of the detector was applied to 
an operational amplifier (PAR Model 215) in the current 
amplifier mode and ac coupled input. The operational 
amplifier output was connected to the ac coupled PAR 225 
low noise preamplifier before going to the ADC. The 
bandpass of the preamplifier was set for 300 Hz to 10 kHz. 

Photomultiplier tubes (1P28 and R166 solar blind) 
were used for visible and UV measurements, The standard 
high voltage dynode chain circuit for photomultiplier 
tubes is for DC applications (Figure 38(a)). This was 
modified (69) for operation of the PMT for signals of 
high frequencies (Figure 38(b)). The PMT circuit was 
normally operated with a negative high voltage of 500V. 
Dl to D9 represent the nine dynodes of the PMT. The 
anode current drains through Ry to ground providing a 
voltage output. This was applied to the ac coupled low 
noise preamplifier, then to the PAR 215 operational 
amplifier in the voltage amplifier mode with gain of 10 
and finally to the ADC. The response of the PMT 


circuit is down 3dB at f=1/(20R,C; ). Cy is the coaxial 
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Figure 38, Photomultiplier tube schematic circuit for 
(a) DC operation, (b) AC operation. 
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cable capacitance and is about 100 pf. For response to 
around 30 kHz, R, was chosen to be 50k ohms. 

5) Sampling Rates 

It is necessary with a versatile Fourier 

transform spectrometer system to be able to acquire data 
at various sampling rates, The importance of sampling 
rate and its relation to the frequency axis has been 
discussed in Chapter III. Most of the measurements are 
now taken using the direct sampling rate of the He-Ne 
laser fringes. Where the number of points for the Fourier 
transform is limited by the size of the computer memory 
(as in the PDP-8/e system used in the early development 
of this work and described in the following section) it 
is desirable to use a coarser sampling interval to 
optimize spectral coverage and resolution. These coarser 
sampling intervals ($2 and $4 sampling rates) are achieved 
by software. Data is taken into the computer at the 
direct rate of the laser fringes. If the #2 rate is 
desired, the program is instructed to signal average every 
other point or for the $4 rate, every fourth point. 
This, however, does mean that 2000 data points must be 
acquired in order to time average a 500 point inter- 
ferogram at the #4 rate. 

The question arises as to why a hardware approach, 
simply a couple of flip-flops on the laser digital clock 


to divide by two or four, could not be used, Time 
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averaging could proceed directly and with simpler software, 
However, this is not possible in our system for a subtle 
reason. The ADC clock must be phase locked to the white 
light and signal interferograms. The laser digital clock 
satisfies this constraint. Figure 37(b) shows 20 
superimposed repetitive scans of the laser digital and 
white light analog signals. The laser digital maintains 
the same phase with respect to the white light signal. 
But any clock signal generated from it by a modulo-2 or 
modulo-4 counter cannot be phase locked to the white 
light and signal interferogram because the state of the 
counter cannot be reproducibly set at the start of data 
acquisition (the start pulse is imprecise with respect to: 
the white light interferogram). 

It is also possible to take data at sampling 
intervals finer than that provided by the direct rate of 
the laser digital clock. This is useful for clarifying 
highly aliased spectra as shown in Chapter III, It 
turns out, in contrast to the above situation, that a 
phase locked clock of higher frequency can be generated 
by a hardware approach. 

A phase locked loop can be used for precise 
frequency multiplication (54). It was mentioned in an 
earlier section that the internal loop of the 565 PLL 
between the VCO and the phase comparator was closed by 


connecting pin 4 and pin 5. If this connection is opened 
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and a modulo-K counter (K=2 or 4) is inserted, the loop 
will operate so as to drive the VCO at K times the 
frequency of the input laser analog signal. The VCO 
free running frequency must be adjusted to about twice 
the laser analog frequency before the PLL can lock onto 
it. Once lock is achieved, the VCO output at pin 4 will 
be precisely K times the laser analog signal frequency 
and phase locked to it. This approach has been used to 
generate X2 and X4 sampling rates, 

6) Electronics with the PDP-8/e Computer 

The computer system used for the early 

development of this work was a DEC (Digital Equipment 
Corporation) PDP-8/e. There are only two digital control 
lines into the computer (the other digital inputs are 
used for the data). One is used for the clock ADC flag 
to tell the computer that the conversion is done and the 
digital data is available to the computer. The other 
control line must then be used to start and also to stop 
taking data, 

The interferometer electronics used with the PDP-8 
is shown in Figure 39 and the timing sequence in Figure 40, 
The circuit for the laser channel and phase locked loop 
feedback control is basically the same as described earlier. 
The start flag, in this case, is a very short pulse 
(about 50 usec) from MONO-2 generated from MONO=-1 and the 


sync. pulse. The white light fringe triggers MONO=3 
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Schematic diagram of the control signals’ 
timing sequence for the circuit in Figure 39, 
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» 


wa 


which opens the modulo-N counter. The overflow pulse 
triggers MONO-4 which gives another very short pulse 
(about 50 usec) which is the stop flag. The counter 
reset used here is more crude than in the latest circuit 
(Figure 34). The counter gate, MONO-3, is adjusted to 
close in a time period longer EAST counts but shorter 
than 2N counts to prevent multiple stop flags. Since 
only one control line is open for the start and stop 


flags, they are combined by a multiplex circuit (OR gate) 


to form one input into the computer. The computer program 


is instructed to sense the first pulse as a start flag 


and the next pulse as a stop flag. 
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D. Computer Processing 


1) The Computer Systems 


The first computer system available was a 
PDP-8/e from DEC, with 16k of core and a DEC-tape based 
OS/8 operating system (Figure 41). The associated 
digitizing system consists of a sample and hold module 
with a 50 nsec aperture time coupled to a 10-bit ADC 
with a 10 usec conversion time, Assembly language 
(SABR) commands were used to input data from the ADC and 
output data to the digital-to-analog converter via a 
DR8-EA 12 channel buffered digital I/O (input/output) 
board. All other programming was written in Fortran 
including the floating point FFT. The maximum number 
of points that can be transformed with this system 
sr O12. 

The present minicomputer data acquisition and 
processing system used with the interferometer is based 
on a DEC PDP=11/10 with 32k of memory and an LA-36 
Decwriter terminal (Figure 42). The RT-11 (real time) 
operating system uses an RK-05 cartridge disk as the 
main storage medium, Accessory storage is provided by 
a TA-1l1 dual cassette unit. Versatile program editing 
and graphics capability is provided by a VI-1ll graphic 
display processor with a VR-17 CRT. Hardcopy output 


is provided by a Zeta Research Series 100 incremental 
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Figure 41, Minicomputer data acquisition and 


processing system based on the PDP-8/e. 
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Figure 42, Minicomputer data acquisition and processing 
system based on the PDP-11/10. 
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plotter (Model 130-5). 

The Michelson interferometer is interfaced to the 
computer by the DEC laboratory peripheral system (LPS) 
which contains the ADC and the sample and hold amplifier. 
The analog signal interferogram is connected to one of 
the eight channels of the LPS multiplexer and preamplifier 
before going to the sample and hold and 12-bit ADC 
(20 usec conversion time). The laser digital clock, the 
start and the stop flags are connected to the LPS 16-bit 
digital I/O, 

The data acquisition program was written in 
assembly language (MACRO). All other programming was 
written in Fortran including the floating point FFT, 
graphics and plotting programs. The number of points 
that can be transformed is 4096, All the spectra shown 
in this work were calculated from 4096 point inter- 
ferograms, unless otherwise stated. The interferograms 
were fully double-sided excent those shown in Chapter V. 
Also the spectra were calculated as the amplitude spectra 
without using any phase correction techniques (except 


the spectra in Chapter V). 
2) Computer Programming 
The programs used to process the inter- 
ferogram and manipulate the data are presented here in 
descriptive form. The actual programs and detailed 


commentary are presented in the Appendix. Figure 43 
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Program sequence for data 
acquisition and processing. 
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gives the program flow to process the data, 

DATACQ is the data acquisition program. The actual 
acquisition subroutine is written in assembly language 
(MACRO). The program begins to acquire data upon arrival 
of the start flag with one point being stored for each 
laser digital clock pulse, The stop flag halts data 
acquisition. The raw interferogram is then summed 
(sequenced from the last data point) into an array which 
is double word to avoid addition overflow. The program 
then waits for the next start flag and the sequence repeats 
until the required number of scans have been taken. Then 
the summing array is divided by the number of scans and 
converted to single word to give the time averaged inter- 
ferogram. The interferogram data points are written onto 
the disk as the file FIN1 in integer format. 

PLTIN1 converts FIN] to a floating point file 
which is stored as FTIN2, FINI] and FIN2 contain 
additional information about the parameters used, such 
as the number of points and scans, date and title, and 
the maximum and minimum values in the interferogram and 
their point locations. 

APODL performs the apodization of the interferogram 
stored in FTIN2, The interferogram is multiplied by a 
mathematical function which can vary from zero to one. 
There is a choice of functions: boxcar, triangular, and 
Gaussian. Others could easily be programmed. The degree 


of apodization is determined by the choice of the indices 
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in the case of boxcar and triangular functions. For the 
Gaussian function, the parameter FAC is entered to 
determine the width of the function and Z to determine 
ae center point of the function (usually 2048 for a 
4096 point double-sided interferogram), The apodized 
interferogram is stored in the file FTINLO, 

COOLF4 performs the fast Fourier transform on the 
data in FIN1O and outputs the unscaled amplitude spectrum 
onto the disk as the file FTN12, 

Both the interferogram and amplitude spectrum can 
be viewed on the CRT with the programs INVT and SPVT 
respectively. Any portion of these data can be displayed 
and expanded, Listing of the data points on the terminal 
can also be obtained. 

Hardcopy can be obtained on the incremental plotter 
of the interferogram, the apodized interferogram, and 
the spectrum with the programs PLTIN2, PLTINA, and TEMP 
respectively. Any portion of the data can be plotted 
with any desired expansion. The spectrum can also be 
scaled to give expanded amplitudes of small peaks with 
the program PLOTSP. 

PKLO searches for peaks in the spectrum above a 
specified threshold value and outputs their point 
location, wavenumber and wavelength (for the specified 
aliasing region). The program also calculates the 


scaling factor so that absolute peak amplitudes can be 


r - e : ma i i 
, 7 4 7 y an 7 - ‘ ‘ i 
i apes Toy iin Vie a) Agee . es 
gee ng a ’ <hAcy Oils) ) ee te on 4 
“7 : ven! ( 
f ’ ) rs 
A! : Ww m 7 i A ' 7 ~ A» Py A 
Pettey, Mge> peeeto 
| 
ne 
> i, d 
| Ww): Qt te ie cies 
: t 
Tt 4 ~ ¢ oi 
y re 4 4 j } ; 
f 7 f 7 
. te if hey + 4 ly 
I 7 ios 


: rt 
: i a tos ncn 
; ‘ 
“Viel pene a ol 
. 7 mam 7 
t - = a 
—— 
' g nas 
! 
hod is 
= He 
» 
ho 
‘ 
y > I 
\ 
i ; ¢ “sft 
=r - 
4 M4 " ry a 
} ay i — ve 
Bi) . ~ , 
. 
i 
P 
1 i 
1 | 
‘ A ag il 4 ws 
i eee fae 4 
. ve f bg ey 
Cre 
= ‘ol : ed 
j ¢ ‘ AP win” iY “3 Le ef 
S 7 
Thi: vy 
Are 
ce \ oo 
X bd MD “a 
™, : 
5 -_ 
im “1h oo és 
oat wt 74 ~~ “= 
: r 
P foal A 
H ! . a +p! 34 a 4 . 


1 
Y 
raed oe: 
eS 7 7 
at 7 * oe 


22 


obtained for quantitative measurements, 

One place where further program development is 
needed is in the number of points the Fourier transform 
program is capable of processing. At present, the computer 
memory of 32k on the PDP-11/10 limits the "in-core" 
transform size to 4096 points (floating point FORTRAN). 
To improve spectral resolution, longer interferograms 
(larger optical retardation) must be transformed. There 
are procedures for performing the Fourier transform in 
segments using the disk memory for temporary storage, 
therefore, allowing much longer interferograms to be 
processed by this computer. Another possibility to 
achieve this is to have a direct communications line from 
the PDP-11/10 to a large central computer, The mini- 
computer would perform the data acquisition, then, the 
interferogram would be transmitted to the large computer 


which performs the Fourier transform. 
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CO ACDERy Vila 
SPECTRAL MEASUREMENTS WITH THE 
FOURIER TRANSFORM SPECTROMETER 


The Michelson interferometer system was described 
in detail in Chapter VI. Earlier chapters discussed 
special factors that had to be considered in extending 
the use of Fourier transform spectroscopy to atomic 
spectrochemical measurements, A few spectra were shown 
to illustrate important points in the earlier chapters. 

In order to explore the capabilities of the 
Fourier transform spectrometer that was built, various 
spectral sources in different spectral regions were 
measured, These include measurements from the near-IR, 
visible, ultraviolet, and mid-IR spectral regions. 
Representative spectra are shown to point out the 


versatility and capability of the spectrometer. 
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A. Near Infrared Measurements 


The flame emission spectra of the alkali metals, 
Li, Na, K, Rb, and Cs, fall in the 600 to 900nm region, 
One detector which has response in this spectral region 
is the Si photodiode. Its maximum sensitivity is around 
lum, The detector is mounted in an aluminum housing 
which fits into the detector mount of the interferometer, 

The instrumental setup for the flame source 
consists of a Varian Techtron manual gas control box, a 
Varian nebulizer and a premixed burner. The slot burner 
head was removed and a cylindrical flame head (70) was 
constructed, The burner head consists of an array of 
stainless steel capillary tubes of 1 mm inside diameter. 
The capillary tubes are stacked side by side and held 
together by epoxy inside the outer cylinder. The top 
surface of the array of about 60 capillaries is flat and 
the array is about 11 mm in total diameter, 

The flame is then actually a composite of small 
flames, one from each capillary tube. The blue inner 
cones of the flame are 1 to 3 mm high, but the rest of 
the flame extends about 30 cm in height. The air/ 
acetylene flame is very stable. For the alkali metals 
a cool flame is sufficient for excitation. The manual 
gas box was set for flows of air=7 units and 


acetylene=2.5 units, This corresponds to 12,3 and beta! 
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liters/minute respectively. The burner was placed at the 
focal point of the input off-axis parabolic mirror to 
collimate the light to the interferometer. The obser- 
vation height for the atomic emission was set by 
rotating the mirror so that the focal point was about 

5 cm above the top of the burner, 

Figure 44 shows a spectrum of the five alkali 
metals listed above. The spectrum was obtained by taking 
the Fourier transform of the interferogram resulting 
from a solution of the five elements at about 
100 parts per million (ppm). The solutions were made 
from the chloride salts of the elements, Fifty consecutive 
scans were averaged and the interferogram was 
apodized with a Gaussian function (FAC=2.5). The 
sampling rate used was the direct rate derived from the 
He-Ne laser. The digitized interferogram is shown in 
Figure 45 where the middle 3000 points of the 4096 point 
double-sided interferogram is plotted. To illustrate 
again the point made in Chapter V about phase errors, 
note that the interferogram is quite asymmetric. Yet, 
the spectrum is not distorted because it is calculated 
from a double-sided interferogram, The Michelson 
moving mirror drive was set at about 0.2 Hz or 5 seconds 
per scan, The 4096 points of data on one scan were 
acquired in a 1 second period during the 5 second 


mirror drive cycle. Therefore, although about 4 minutes 
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Digitized interferogram which produced the alkali metals 


spectrum in Figure 44, 


Figure 45, 
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are needed to acquire the data for 50 scans, the actual 
measurement time is only about 1 minute, 

It should be noted that, since the sampling 
interval was 0,6328 um, all the lines in Figure 44 are 
aliased as explained in Chapter III, The wavenumber 
axis is labelled for region (5-8) (see Table II). 

The Li, K, Rb, and Cs peaks lie in this region, but the 
Na doublet is actually aliased in from region (9-12). 

The measured wavelengths are tabulated and compared with 
the literature values (71) in Table IV. The agreement 

is quite good. It should be emphasized that a spectral 
range of 300 nm is covered in this spectrum, illustrating 
the wide simultaneous spectral coverage possible with 
Fourier spectroscopy. 

Quantitative aspects were studied with the lithium 
peak at 670.78 nm. An analytical curve for lithium 
flame emission is shown in Figure 46. The lower 
concentration limit is about 0.1 ppm at which point the 
S/N for the lithium peak decreases to about three. The 
slope of the straight portion of the log-log plot 
ranges from 0.95 to 1,05 for all the analytical curves 
that have been run. The slope in Figure 46 is 1.05. 
This shows that the Fourier transform spectrometer is 
capable of quantitative measurements. 


The spectra for two of the points on the analytical 


curve are shown in Figures 47 and 48 for 5 ppm and 0.5 ppm 
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TABLE IV 
Wavelengths Calculated From the Alkali Metals 


Spectrum Shown in Figure 44 


aaa Ee a OO 


Element Literature (mm) Experimental (nm) Anm 


Li 670.78 670.63 Ogs 
Na DOI. 09 589.47 Oe 
589.00 588,87 Oa 

K 766.49 766432 On l7 
769.90 769.74 Oe6 

Rb 780,02 IQS 0.09 
794,76 794,64 0.12 

Cs Spe. Soleo? 0,12 
894, 35 894,17 0.18 


} in 2 ) 
em Aion cee ale gn ii eth woe 


i | 
~— 


4 ben Lo) hae 
Ee. #08 


500 


100 


UNITS 


ARBITRARY 


Figure 46, 


Jc eek Were lt EL er Tees TT el Lee ROSE EMRE Beeeel Uae tag ter ps (Aa eit See mea paral e. J BS 


1 10 100 
LITHIUM CONC. (PPM) 


Analytical curve for lithium flame emmission. 


, mate , a 
a 
_ Wwevwrar ‘ ee ee eee 
v1 oN ; 
hy ;/ 


rat) OMI) se) IWETNE 


| r . ; 
‘owe eal ¥ wirisfall au evita pa ‘VPA le 


131 


71062 


82086 


*uUOTINTOS umTyaAtTT wdd ¢ e Fo umaqoeds 


YSSNNN SAV 
gZ0EIl 8-Z08EL 


mu g/°0L9 FI 


°/y eANnsty 


8208S! 


—— + - 


AGNLNdNV 


X a joa a a a TS were Tere 


| See ia 


1 


j 
ro) én ® 
7 Tn 7 


Ty oar) 


Sy 


7062 


*uoTANTOS unTYyaTT wdd ¢*og e FO wumaqoeds 


YASWNNIAVM 
8° Z08IL 8 2O8El 


Fada Ra PY 


wu g/°0L9 


°gy eansTy 


FT 


EA RE OS, 


: aaa It 
S25 ep caer ae 9 Joe a a ie ean 


87208Sl 


FJANLIIdNV 


lithium respectively. The interferograms were double- 
sided and were each the average of 50 scans. Gaussian 
apodization was applied to the 4096 point interferograms 
(0.6328 um sampling interval). The absolute peak 
amplitude was calculated from the scaling factor needed 
to scale the spectrum from 0 to 4096 on the vertical 
axis. The band spectra around 10800 to 9800 em? > are due 
to water bands around 960 nm (72). The amplitudes of 
these bands remain relatively constant. The bands in the 
5 ppm Li spectrum. appear smaller because both spectra 
have been scaled between 0 and 4096, but, the Li peak is 
ten times larger for the 5 ppm spectrum than for the 0.5 
ppm spectrum. 

To extend the measurement of spectra to shorter 
wavelengths, studies were carried out in the visible 
region. The atomic spectra of most elements lie in the 


visible region of the spectrum. 
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B, Measurements in the Visible Spectral Region 


Measurements made in the visible spectral region 
used a photomultiplier tube (1P28) as the detector. The 
negative high DC voltage for the PMT was supplied by a 
Heath EU-/701-30 power supply module and was set for 400 
to 500 V. The dynode chain for the PMT was mounted inside 
the PMT housing. The housing has an adaptor piece which 
fits into the detector mount of the interferometer, 
Different size apertures could be placed on the adaptor 
Piece. Most of the spectra were taken with a 1.0 mm 
aperture, 

Sources measured were flame emission of Ca, Sr, Ba, 
Na, Li and hollow cathode lamps (HCL) of several elements. 
All the spectra are from double-sided interferograms 
apodized with a Gaussian function (PAG=245) seine drrect 
sampling rate (0.6328 um) was used and 4096 points were 
acquired. 

Hutiyascansswere averaged, This number was 
chosen as a compromise between improved S/N ratio with 
more scans and shorter time of data acquisition with 
fewer scans, To test the S/N ratio improvement with 
more scans, the baseline noise of Na spectra was 
measured for interferograms of l, 10, and 100 scans. 

The increase in the S/N ratio (decrease in noise with 


the signal constant) roughly follows the square root 
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of the increase in the number of scans, 
1) Atomic Flame Emission 

The flame emission spectra of Lay, Na, sBa 
Sr, and Ca is shown in Figure 49, The concentrations 
emer ba SOO °ppm, Sr 100, Na 1, Li OF anda ceulOor 
The flame used was the same air/acetylene flame used for 
near-IR measurements, No attempt was made to optimize the 
flame conditions, Measurements using a nitrous oxide/ 
acetylene flame would give better emission signals for 
many elements including Ca, Ba and Sr. Another excellent 
source for emission is the inductively coupled plasma. 
This is a very good source for simultaneous multielement 
analysis, 

The barium peak is overlapped by a Ca0OH 
band centered at about 555 nm. Again, the spectral lines 
are all aliased. The labelled axis is for region (9-12) 
which includes the lines of Na, Ba, Sr, and Ca. The Li 
line is in region (5-8). The interferogram is shown in 
Figure 50 (the middle 3000 points are plotted). Note 
that since the Na doublet dominates the spectrum in terms 
of peak amplitudes, the interferogram has the general 
shape of the Na interferogram (a beat pattern from the 
summation of two cosine waves of slightly differing 
frequencies). The interferogram is fairly symmetrical. 


The symmetry is very dependent on the proper alignment 


of the Michelson mirror. The literature and measured 
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tized interferogram which produced the 


spectrum in Figure 49, 
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Figure 50. 
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wavelengths are tabulated in Table V for comparison. 
Again, the agreement is quite good. 
2) Hollow Cathode Lamp Spectra 

The ability to measure a spectrum with a 
number of spectral lines was studied using hollow cathode 
lamps as the source, A Heath EUW-15 power supply was 
adjusted to give 20 to 35 ma of current (depending on 
the lamp used) to the hollow cathode lamps. The lamps 
were clamped in place such that the cathode was at the 
focal point of the input off-axis parabolic mirror. All 
the interferograms were averages of 50 scans except the 
multielement steel HCL interferogram which was 250 scans. 
Actually, 50 scans of the steel HCL produced good spectra, 
but 250 consecutive scans were acquired to test the signal 
averaging capability of the system (that is, the stability 
of the precise stop flag). 

The spectrum of a multielement steel HCL was 
measured (Figure 51). The elements present in the HCL 
are: iron, chromium, nickel, copper, cobalt, and 
manganese. An optical filter (Corning filter number 7-59) 
with a peak %T at 370 nm and a greater than 10ZT from 
305 to 465 nm was used to limit the spectral bandwidth 
to approximately one aliasing region. Actually, the 
filter used allowed a portion of another aliasing region 
to be measured. The chromium triplet (427 nm) shown at 


the far right of the spectrum in Figure 51 belongs in 
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TABLE V 
Wavelengths Calculated From the Spectrum 


Shown in Figure 49 


LS SS ess SSS SSE 


Element Literature (nm) Experimental (nm) Anm 
Li 6707.78 670,54 0,24 | 
Na 309.29 589.40 O5bo 

589,00 588.80 0.20 
Ba DOI. D2 553,34 0,21 
or 460.73 460.56 Oly 


Ca 422.67 422.49 0.18 
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region (9-12). The wavenumber axis shown on the spectrum 
is for region (13-16) which contains the remainder of the 
spectral lines, 

The entire 4096 point spectrum is plotted for the 
steel HCL to emphasize the density of spectral information 
that is present, An important point to realize is that 
the entire spectrum was acquired in about 25 minutes 
including the Fourier transformation. If 50 scans were 
acquired instead of the 250 scans used for this spectrum, 
the time involved would only be about 6 minutes. 

As a comparison, a scanning monochromator was used 
to measure the steel HCL spectrum, The scanning rate 
was set to give approximately the same resolution as in 
the spectrum taken with the Fourier transform spectrometer, 
The portion of the spectrum using the monochromator 
around the manganese triplet at 403 nm is shown in 
Figure 52, The measurement system of the Heath EU-700 
Series Instrumentation for Spectroscopy was used. It 
consists of: EU-703-62 hollow cathode lamp power supply, 
a 0.35 meter monochromator, EU-701-30 PMT module, 
EU-20-28 log/linear current module, and a EUW-20A 
servo-recorder, The entrance slit width was set at 
50 um, and the scan rate of the monochromator was set at 
0.1 A/sec. Therefore, to measure just this 1.2 nm 
portion of the spectrum required 2 minutes. To cover 


the 305 to 465 nm range in Figure 51, the scanning 
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Spectrum of Mn 403 nm triplet using a 


Figure 52, 
scanning monochromator system. 
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monochromator would require about 260 minutes, The Heath 
system has a rather slow response recorder. The S/N ratio 
was not equalized with that in Figure 51, Even with a 
faster response system, the monochromator is still over 
an order of magnitude slower in acquiring the data. This 
comparison is useful to appreciate the resolution that 
can be achieved along with good speed of acquiring the 
data with the Fourier transform spectrometer, 

An expanded portion of the steel spectrum measured 
with the Fourier spectrometer is plotted in Figure 53 
showing the Mn triplet at 403 nm, The wavenumber axis 
is labelled for the Mn. The chromium triplet at 427 nm 
belongs to the next aliasing region. Expanded plots of 
several other portions of the steel spectrum in Figure 5l 
are shown in Figures 54-56. Figure 54 shows the region 
around the iron 371.99 nm line, Figure 55 shows the 
chromium triplet at 360 nm and Figure 56 shows the region 
around the nickel 352.45 nm line. These expanded plots 
show the large amount of information that is present in 
the entire spectrum, 

Several spectral lines in Figure 51 are tabulated 
in Table VI. The agreement with the literature values 
is very good, Only a few representative lines are listed. 
In excess of fifty lines have been identified with the 
wavelength agreements in the same range as those listed 


in the table. Identification of the lines is 
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Figure 53, Expanded portion from the steel spectrum in 
Figure 51 showing the Mn 403 nm triplet. 
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Figure 54, Expanded portion from the steel spectrum 
in Bea Sl in: thé; Vicinity. of the 
Fe 372 nm line. 
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Figure 55, Expanded portion from the steel spectrum 
in Figure 51 in the vicinity of the 
Cr 360 nm triplet. 
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Figure 56, Expanded portion from the steel spectrum 
in Aa has 51 in the vicinity of the 
Ni 352 nm line. 
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Wavelengths Calculated From the Spectrum of the 


Element 


Literature (nm) 


403,449 
403,307 
403.076 
374.550 
371,994 
360.532 
359.348 
357.868 
352.454 
5216 390 
324,754 


TABLE VI 


Experimental (nm) 


403,443 
403.286 
403.066 
374,548 
371.968 
360.505 
559,590 
357.842 
352.440 
SZ ioe 
324,726 


Multielement Steel HCL Shown in Figure 51 


a 


4nm 


0.006 
0,021 
0.010 
0.002 
0,026 
0,027 
0.017 
0.026 
0.014 
0.024 
0.028 
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straightforward because of the accurate wavelengths. 
With the Fourier transform spectrometer, the wavenumber 
axis is inherently accurate (depending only on the 
accuracy of the He-Ne laser and the accuracy of the 
sampling). There is no external calibration necessary. 
The spectrum of other hollow cathode lamps have 
also been measured to further illustrate the capability 
of the Fourier transform spectrometer, Figure 57 shows 
the spectrum of a calcium, aluminum, magnesium three 
element HCL, The last 3000 points of the 4096 point 
spectrum are plotted in Figure 57(a). Again, the Corning 
Filter 7-59 was used. The wavenumber axis is for region 
(13-16). The calcium 422.67 nm line belongs in another 
aliased region (9-12). Figure 57(b) shows an expanded 
portion of the spectrum, The Mg lines around 383 nm, the 
Ca line at 393.37 nm, and the aluminum lines around 
396 nm are shown. The smaller peaks in the spectrum are 
due to the neon filler gas. Table VII tabulates the 
Calculated lines from four independently measured spectra. 
The spectrum of a chromium HCL is shown in Figure 38, 
The optical filter (7-59) was used. Only about 3000 
points of the spectrum are plotted in Figure 58(a). The 
wavenumber axis is for region (13-16). The Cr triplet 
at 360 nm is expanded in Figure 58(b). The Cr 427 nm 
Eciplet at thestarerighteingragure 58(a) is expanded in 


Figure 58(c). This triplet belongs in region (9-12). 
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Figure 57.a)Spectrum of a Ca-Al-Mg hollow cathode lamp. 
b)Expanded portion of the spectrum in a) showing 
= some magnesium, calcium, and aluminum lines. 
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Spectrum of a chromium hollow cathode lamp. 
Expanded portion around the Cr 360 nm triplet. 


(c) Expanded portion around the Cr 427 nm triplet. 
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TABLE VIII 
Wavelengths Calculated From the Spectrum of the 


Chromium HCL Shown in Figure 58 


Element Literature (nm) Experimental (nm) anm 


Cr 357,87 Boys OF 0.00 
Cr 659935 BD IG30 0.01 
Cr 360.53 360.53 0,00 
Cr 428,97 428,92 0.05 
Cr 427.48 427.43 0,05 
Cr 425.43 425.40 0.03 
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The axis is relabelled for this region in Figure 58(c). 
The measured wavelengths are compared with the literature 
values in Table VIII. The agreement is excellent, 

3) Effect of Large Spectral Features on Small Peaks 

7 One disadvantage of the Fourier spectroscopy 

technique is that, if one is interested in a small spectral 
peak, large unwanted signals may have to be measured as 
well, Optical filtering may help to reduce the unwanted 
signal, but there may be instances where the two signals 
are too close in wavelength to allow effective optical 
filtering. As the unwanted signal increases in amplitude, 
there may be a decrease in the S/N ratio of the small 
peak of interest, This effect was investigated by 
measuring the flame emission of a strontium/calcium 
solution and a strontium/sodium solution. 

The flame emission of a solution of 10 ppm Sr and 
1000 ppm Ca was measured (Figure 59). The Sr 460.73 nm 
line occurs near the 21802,.8 om”! mark, The Ca 422.67 nm 
line at the right and the CaOH bands to the left are much 
larger than the Sr peak, yet, the Sr signal is still 
measurable, Figure 60 is an expansion of the Sr peak. 
Note the baseline noise in this and subsequent expanded 
plots of the Sr peak in this section. The noise levels 
will be discussed shortly. 


Another example is shown in Figures 61 to 65 where 


the Sr concentration was kept constant at 10 ppm and the 
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Figure 61, 
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Flame emission spectrum of a solution of 
10 ppm strontium by itself. 
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Expanded portion of the Sr 10/Na 10 ppm 
spectrum showing the Sr peak. 
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Na concentration was increased, In Figure 61, the spectrum 


was obtained with Sr by itself, Only the region around 
the Sr peak is plotted. Figure 62 shows the spectrum of 
Sr at 10 ppm and Na at 10 ppm (only the first 3370 points 
are plotted, the remainder of the points being continued 
baseline). The Sr peak is at the right and is fairly 
small, Figure 63 is an expansion of the region around 
the Sr peak. Note that the noise level has increased 
very slightly from Figure 61 to Figure 63, Figure 64 
shows the spectrum of Sr at 10 ppm and Na at 100 ppm. 
The Sr peak can barely be seen. (The peaks around the 
Na doublet are images due to reflections off the beam- 
splitter substrate and compensator surfaces and are 
dependent on the mirror alignment.) The portion of 

the spectrum around the Sr peak is expanded in Figure 65. 
The noise level now has increased considerably. 

These two examples, however, lead to differing 
views of the effect of large spectral features on small 
peaks. The Na/Sr example shows a decrease in the S/N 
as the large Na peaks are increased in amplitude (10 ppm 
Na to 100 ppm Na). This is as expected when using a 
signal noise limited detector such as a PMT, where the 
noise increases as the square root of the signal intensity 
(shot noise). Therefore, if much unwanted signal is 
present, the baseline noise increases ( in Fourier 


spectroscopy, noise in the interferogram translates into 
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a noise that is distributed throughout the Spectrum). 

But, the Ca/Sr example seems to contradict the 
above view. From inspection of Figure 59, the Ca0H 
bands are broad and intense, therefore, there must be 
orders of magnitude more total signal than in the Na/Sr 
spectra, Yet, the Sr peak is not orders of magnitude 
noisier than the Sr peak with just Sr present (Figure 61), 

The above contradiction cannot be explained without 
further detailed measurements of signal and noise. There 
must be other factors which contribute to the baseline 
noise, One point to note is that, in addition to the 
varying flux of light as optical retardation is changed 
(this is the interferogram), there is a constant flux 
present as well, This may contribute to the noise since 
as the overall intensity of incoming light increases, this 
constant flux level also increases, The detector receives 
a higher overall flux, hence a greater detector noise 
results. Another factor is that the amplitude spectrum 
is calculated, thus, the baseline noise is always positive. 
Hence, small peaks may get buried in the positive base- 
line noise, Quantizing noise was considered. But 
measurements made with the same preamp gain did not change 
the results to any great extent. 

The behavior of the S/N ratio in the two examples 
is still a mystery. The explanation may be in one of the 


above factors or in others, therefore, further studies 
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must be undertaken, Whatever the outcome and in spite of 


the noise contributions above, the ability to measure 
small peaks in the presence of large spectral features 
is still quite acceptable. The small Sr peak that is 
expanded in Figure 65 is only 0.3% of the Na peak 
amplitude in Figure 64, 
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C, Measurements in the Ultraviolet Region 


The measurements made in the ultraviolet region 
used a solar blind photomultiplier tube (R-166). The 
spectral response of the PMT limits the spectral region 
that is measured to the aliasing region (17-20). A 
magnesium hollow cathode lamp was used as the source, 

Figure 66(a) shows the spectrum obtained from a 
512 point double-sided interferogram using a 0.6328 um 
sampling interval. The resolution is not good enough to 
resolve several of the peaks, The peak numbers refer 
to the lines listed in Table IX. Peak 3 is actually three 
lines and peak 4 is a quintet. To achieve higher 
resolution, the sampling interval was increased to 2,5312 
um (still 512 data points). The lines are now highly 
aliased (Figure 66(b)). The wavenumber axis is shown for 
peaks 2, 3,.and 4 in region 19 and for the lange 
285.21 nm peak which is in region 18. 

The spectral region in Figure 66(b) around the 
peaks 1 to 4 was interpolated using the zero-filling 
technique (73). The result is shown in Figure 6/7. All 
the peaks are resolved and their wavelengths are listed 
in Table IX, The agreement with the literature values 
is very good. 


When the new PDP-11/10 system was available, the 


Mg spectrum was measured again. The interferogram 
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(a) 
Mg 285.2} 
3 
{| 4 
5 
Fe SS 
31,605 em”! S516 Sotch a 391507 ema! 
(316.4 nm) (281.2 nm) (253.1 nm) 
b 
(b) , 
Mg 285.21 
3 
2 4 


CG Se 


=| 
35,556 cm! Region 19 (Peaks 2,3,4&5) 37,531 cm 
35,556 cm"! Region 18 ( Paak 1) 33,581 cm"! 


; i & f a Me hollow cathode lamp (512 point 
eee BEN cay with : 0.6326 um sampling interval, 
(b) with a 2.5312 um sampling interval. 
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TABLE IX 


Shown in Figure 67 


Experimental (nm) 
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Wavelengths of Magnesium Spectral Lines 
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is now 4096 points long. The resultant spectrum using 
the 0,6328 um sampling interval is shown in Figure 68. 
Only the section around peaks 2 to 4 is plotted. This 
spectrum can be compared with Figure 66(a) to see the 
increased resolution with the longer interferogram. 
All the peaks are resolved in Figure 68, The 285 nm 


line is past the left end of the plotted spectrum, 
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D. Infrared Measurements 
Seetclia oliglir lise Ald Scales sce er) 


The Fourier transform spectrometer is also capable 
of infrared measurements. The modularity of the spec- 
trometer allows the beamsplitter assembly, the source, 
and the detector to be easily changed for mid-IR 
measurements, The beamsplitter assembly uses NaCl 
substrate and compensator instead of quartz. Infrared 
studies were not done in detail. Measurements were 
taken to show the.IR capability only. 

A globar (air cooled) was used as the source. A 
polystyrene absorption spectrum was obtained by placing 
a film sample directly in front of the TGS detector in 
the detector assembly. 

The background spectrum is shown in Figure 69 (a) 
and the polystvrene spectrum in Figure 69(b). Only the 
mid-IR region of the spectrum is plotted. 

The double-sided interferograms are shown in 
Bigure (0 (background) and Figure 71 (polystyrene) to 
illustrate the difference in the interferograms between a 
broadband source and the line emission type sources 
presented earlier in this chapter. The middle 3000 
points of the 4096 point interferograms are blotted. 

Infrared measurements are the subject of a future 


studv and are therefore only mentioned to show that the 


Spectrometer can be used in the mid-IR. 
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Figure 69, (a) Background spectrum with a glowbar 
source in the mid-IR region. 
(b) Polystyrene absorption spectrum. 
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CHAPTER VIIT 
SUMMARY 


The Michelson interferometer described in Chapter 
VI has developed considerably from the initial design. 
The earlier interferometers were slow scanning using a 
micrometer drive. The interferogram was recorded with a 
digital voltmeter, The raw data was then translated 
onto computer cards and taken over to the computer 
center for the Fourier transform calculation and output 
plots of the processed data, The present design is 
fast scanning and is laser fringe referenced for 
digitization. The data is acquired directly and processed 
immediately with the minicomputer. 

The interferometer is now at the stage where 
emphasis can be placed on experimental applications 
rather than on instrumental developments. This does not 
mean that further improvements and changes cannot take 
place. However, the measurements taken to this date 
show that the Fourier transform spectrometer is capable 
of producing medium resolution, high quality spectra. 

To achieve this capability at shorter wavelengths 
than the mid-IR, several factors had to be considered. 


The higher spectral frequencies immediately create diffi- 


culties in the proper sampling of the interferogram. 


. . , t 
Since avoiding aliasing 1s not practical normally, it mus 
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be accepted and even, Perhaps iputsato vse. 

The measurement of atomic line Spectra in the 
UV-visible region differs from that of absorption 
Spectra in the IR. The instrumental line shape and 
the means to change it for narrow emission type signals 
has been studied. The fundamental line shape in 
Fourier transform spectroscopy (sin x/x) introduces 
undesirable features (side lobes) in the spectrum, 
Therefore, apodization of the proper types must be 
applied to the interferogram, 

Another problem with narrow emission type signals 
is the difficulty of phase correction. This was 
investigated in order to find methods of overcoming 
spectral distortions due to phase error. A solution 
was not found to correct for phase problems with 
single-sided interferograms., Therefore, the problem 
was avoided with the use of double-sided interferograms, 

With a better understanding of the factors relating 
to atomic spectra, some measurements were under- 
taken. Sources that were used include flame emission and 
hollow cathode lamps. A wide range of wavelengths was 
measured extending from the ultraviolet to the near-IR 


Spectral regions. Some of the advantages of Fourier 


transform spectroscopy are shown in these measurements, 


These include accurate wavenumber ax1s, wide 


simultaneous spectral range, speed of operation, and 
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good resolution in a compact system. 

The performance and versatility of this Fourier 
transform spectrometer is proven in these measurements 
to be quite good. It can be concluded that the Fourier 
transform technique is an alternative to the standard 
instrumentation for atomic spectrochemical measurements. 
Further developments in automation, price-performance, 
and instrumental improvements are needed to make the 
technique part of the repertoire of standard 


instrumentation. 
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APPENDIX 


The computer programs listed here were used to 


acquire and process the interferograms, Reference should 


be made to Figure 43 and the general description of the 


programs in Chapter VI D. 
The data acquisition subroutine was written 

in PDP-11 assembler language (MACRO). All remaining 
programs were written in FORTRAN IV to run on a 

DEC PDP-11/10 minicomputer. Commentary is sprinkled 
throughout the programs to make interpretation easier. 
A new program listing starts on a new page. The 
program name is given at the beginning of a new 


program listing. 
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NXOTF#1000+/FL0OATONG) 

XO FLOAT ONXOLP 
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PROGRAM NAME "PLYINO* 

INTERFEROGRAM FLOTTING ROUT TAK 

FOR THE ZETA INCREMENTAL. PLOTTER. 

THE FLOATING FOLNT INTERFPEROGRAM TS PEAT 
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biter PLOTTING POUT 
FOR THE ZETA INCREMENTAL PLOTTER 
THE FLOATING POINT INTERFEROGRAM TS 
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CLL. i 


GMs Sy elk 
. NUMA 
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TALL AXIS (0.90. » SHAMIL [TI 


Pic ULSI 


Ree Oe Fede O, L4eDATE oO, Oe 12) 
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AMAXSY C1) 

AMINEY CL) 
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CALL AXES CO. 0.9 CMR Lee 4e CAN AHA ZATNC 9 Oe KOMP TEL KOMP TED) 
(CAL AXIS CO. Oe» AMPLITUDE 7 9 $998. 290.» TONPTHEL) » PENPT HE? DD 
COM Ut KOn Om0s Oy 5) 

ie RTE 

GO YO 4 

NGI 

ENG 


2 
f ‘ 
ara) 
5 
\ 

Lh 
Pad 
: 


{ A) : : Uar'r am | i ; 


: 7 
; ooh 


“> 
A - a 
Ula ae ‘ 
cane im 
Dheeuk we 
ere ane be | 
ope ty sae 
a at wr 
sea 
| toy 1D okt 08 
/ ngs Ab Jaime 
if 
¥ 
th 


Spier ous on es Cone 


Sp ules Bee, 


wO4 


PROGRAM NAME "TEME 

INFUT TO PROGRAM [9 FINL2.naT 
“— ee eee wt Ee pee, pore nen i as oa ie 4 ! 
THE ENTIRE AMPLITIOE ¢ ‘LIM Hi 


9 SCALED 6-494 0 


OUTPUT IS PLOT OF POINTS BETHEE TETED CH—l LIMITS 


THE PLOTTED SECTION TS NOT EXEAN 
IN THE PROGRAM "FLOTSE", THE AMPLTTHnr ean 
0 TO 4094 EXTENDS Ta @ INOMES; ee ee 


ene Ne 
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FORMAT(’ “ENTER CM-1 LIMITS OF ALTASING REGIONS 2F10,2°) 


REANCS SOL) CMl»y CM2 
FORMAT (2F 10,2) 

WPS lie oy TEMS) (AO) In GF 
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Cits-CM2 

CONT TNE 
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NEUE Qoe } 
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REWING 22 
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TF CAMIN «GE. CTD) AMTNeY CT) 
CONTINUE 

S# 4096 .O7 CAMAX AMIN} 
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